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SUMMARY 


A computational method has been developed to provide an analysis for 
complex realistic aircraft configurations at transonic speeds. Wing-fuselage 
conf.gurat.ons with various combinations of pods, pylons, nacelles, and 
wmglets can be analyzed along with simpler shapes such as airfoils, isolated 
wmgs, and isolated bodies. The flexibility required for the treatment of such 
diverse geometries is obtained by using a multiple nested grid approach in the 
fmite difference relaxation scheme. Aircraft components (and their grid sys- 
tems) can be added or removed as required. As a result, the computational 
method can be used in the same manner as a wind tunnel to study high-speed 
aerodynamic interference effects. The multiple grid approach also provides 
h.gh boundary point density/cost ratio. High resolution pressure distributions 
can be obtained. Computed results are correlated with wind tunnel and flight 
data using four different transport configurations. Experimental/computational 
component interference effects are included for cases where data is available. 

The computer code used for these comparisons is described in the appendixes of 
this report. 
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INTRODUCTION 


There has been a dramatic growth in the development of transonic compu- 
tational methods during the past ten year period. This growth is stimulated 
by the important role that high-speed non-linear methods now play in the de- 
tailed aircraft aerodynamic design process. Many advances have been related 
to improvements in computing efficiency. In addition, the ability to treat com- 
plex geometric shapes has steadily improved. Murman and Cole^, Bailey and 
Steger , Keller and South (3) , and Jameson and Caughey (4 > have made import- 
ant contributions to this evolutionary process. The primary goal of these ef- 
forts is the development of tools for solving practical aircraft design and anal- 
ysis problems. 


The first analysis schemes were applied to simple 2-D airfoil shapes and 
axisymmetric bodies (see Fig. 1). Soon after this, isolated wing methods 
evolved. Computations were performed on simple wing-body shapes during the 
mid- 1970’s period. By 1980, aircraft configurations with various combinations 
of nacelles, pylons, and winglets^ could be analyzed. 

Two different approaches have been implemented in the past. The first, 
and simplest, employs a small disturbance- type flow equation coupled with 
mean-surface boundary conditions. In time, it was noted that the classical 
transonic small disturbance equation was not well suited to the analysis of 
swept wing flows. The proper shock strength and the velocity at which the 
flow equation changes type is not inherent in the equation when shock waves 
have appreciable sweep. For this reason, today’s small disturbance formula- 
tions typically incorporate a modified or extended small disturbance flow 
equation. The second, and more sophisticated approach, employs the full po- 
tential equation (or Euler's equations) coupled with surface -conforming bound- 
ary conditions. This approach is complicated, not so much by the use of the 
more complete flow equation, but rather by the difficulties associated with the 
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application of exact boundary conditions. As configurations become more com- 
plex, it becomes difficult to construct a suitable coordinate system. As a re- 
sult, the availability of full potential equation methods has lagged that of the 
small disturbance methods by several years. 

Until recently, small disturbance and full potential equation methods were 
constructed about a single continuous grid system. This approach provides a 
smooth transition from the near to the far field regions when proper stretching 
functions are used, but for a case as simple as a wing-body combination, 
where axial grid point resolution must be high for the body and spanwise res- 
olution must be high for the wing, problems arise. Problems are caused by 
time and core storage restrictions of current computers and inefficiency of the 
combined axial/ spanwise grid system which stretches into the far field where it 
is not required. Transonic wing- body analysis codes with single continuous 
grid systems exhibit poor resolution on either the body or the wing because 
these surfaces are not aligned. These difficulties can be overcome by not 
implementing the traditional single grid approach. 


SPE - SMALL PERTURBATION EQUATION 
FPE - FULL POTENTIAL EQUATION 
EE - EULER EQUATIONS 


1971 - SPE 

1972 - FPE 
1976 - EE 



AIRFOIL 


1976 - FPE 



AXISYMMETRIC BODY 


1972 - SPE 
1974 - FPE 
1981 - EE 



WING 


1975 - SPE 
1977 - FPE 


1978 - SPE 

1979 - FPE 



1980 - SPE 



WING-FUSELAGE-NACELLES-PYLONS- 
WINGLETS-CAN ARDS 


R84-1 1 37-001 D 

Figure 1 Growth in Geometric Complexity for Transonic Flow Simulations 
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Grid orthogonality presents another problem for "exact" boundary condi- 
tion formulations. Boundary/grid orthogonality is easy to achieve for simple 
2-D and 3-D shapes, but for complex shapes such as a wing-nacelle-pylon ar- 
rangement, it may be physically impossible to achieve an orthogonal grid for 
all surfaces simultaneously. This often results in a lack of orthogonality at 
the critical pylon-nacelle station. This could, depending on the severity of 
the flow condition and the geometric shape, produce a jump in the metric 
across the pylon-nacelle span line. The resulting "numerical" interference 
effect might be interpreted as an "aerodynamic" interference effect during 

applications. This also suggests, then, a re-evaluation the conventional single 
grid approach. 

When configurations are examined by wind tunnel testing, various compo- 
nents such as pods, pylons and winglets are removed and reattached to evalu- 
ate aerodynamic interference effects. If a three-dimensional computational 
method is developed using the single grid approach, it will be necessary to 

find a different coordinate mapping for each configuration modification. Per- 
haps a separate computer code would be required for each geometric shape in- 
volved. The configurations boundary points, mesh resolution, and ultimately, 
the computed flow field would vary with each change to the coordinate system. 
This problem is minimized if multiple grid systems are used. A single com- 
puter code can be used for the analysis of isolated bodies, isolated wings, and 
wing-fuselage configurations with nacelles and additional lifting surfaces. 
Boundary surfaces, along with their embedded mesh systems, can be added or 

removed for analysis much in the same manner as wind tunnel model parts dur- 
ing experimentation. 

A simple multiple nested grid approach was described in NASA CR-3243 
(Reference 6). The simplicity of the grid embedding scheme is stressed be- 
cause the probability of successfully simulating a transonic flow about a com- 
plex configuration will be increased if the methodology is simple. This report 
describes the extension of the Reference 6 method to include the simultaneous 
treatment of multiple wing and body surfaces. In particular, the provisions 
that have been made for the analysis of pods, nacelles, pylons and winglets 
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are noted. The small disturbance character of the original wing-body method 
is retained. The mean surface approximation permits the extensive use of rec- 
tangular coordinate arrays. This plays an important role in keeping the meth- 
od simple. 

The basic transonic Wing-Body COde (WIBCO) was described in Reference 
6. As a result of this effort, a new computer code has evolved (WIBCO-PPW) 
which in addition to the wing, body and wing-body capabilities of WIBCO, in- 
cludes the new pod, pylon and winglet analysis capabilities. This report 
should be used along with Reference 6 since certain elements of the two codes 
which are common are not repeated herein (i.e., wing viscous effects, arbi- 
trary fuselage modeling, wing, body, wing-body flow simulations). 

The WIBCO-PPW Code has been developed for both IBM and CDC type 
computers* Comments on usage and features of the NASA CDC version of the 
code have been included as Appendix B. 


* The author would like to thank Dr. Perry Newman for many valuable dis- 
cussions during this code development effort. In addition, Mrs. Kara Haigler 
made significant contributions by refining and verifying the NASA CDC version 

of the code. 
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NOMENCLATURE 


SYMBOLS 

b 

C,c 

C D’ C d 

Cp,Gf 

C L > C l 

SpSn 

C P 

f 

g 

L 

M. A.C. 

M 

MFR 

NPR 

N, n 
P 

R,r 

R 

T 

u,v,w 

V 

X,Y,Z 

a 

6 

r 

Y 

AX, AY, AZ 
A? , An, A 5 
A <f> 


E.n.c 

♦ 


Wing or winglet span 
Wing or winglet chord 
Drag coefficient 
Friction drag coefficient 
Lift coefficient 
Moment coefficient 
Pressure coefficient 
Surface shape function 
Acceleration due to gravity 
Pod or nacelle length 
Mean aerodynamic chord of wing 
Mach number 

Nacelle inlet mass flow ratio 
Nacelle nozzle pressure ratio 
Surface normal 
Pressure 

Nacelle/pod radius 
Gas constant 
Absolute temperature 
Velocity components 
Velocity 

Cartesian coordinate axes 
Angle-of -attack 
Angle-of-yaw 
Circulation 

Specific heat ratio (Cp/Cy) 

Mesh spacing in X, Y, and Z directions 

Mesh spacing in K , n , and z directions 

Change in perturbation potential between iteration cycles 

Boundary layer displacement thickness 

Wing or winglet span position (2 y/b) 

Wing skewed coordinate axes 
Perturbation potential 


Winglet toe-out angle, measured from symmetry plane 
Winglet cant angle, measured from wing plane 


Subscripts : 


av,AV 

Average 

Body,B 

Body 

c 

Computational 

cold 

Cold jet 

D 

Dummy value 

e 

Exhaust 

FIELD 

Field 

FRICTION 

Friction 

g 

Geometric 

hot 

Hot jet 

I 

Induced 

i 

Inlet 

ib 

Inboard 

i,j>k 

Mesh indices 

J 

Jet 

loc, LOC 

Local 

MAX 

Maximum 

0 

Total 

ob 

Outboard 

pod, POD 

Pod 

pyi 

Pylon 

REF 

Reference 

S 

Static 

WAVE 

Wave 

WING,W 

Wing 

x,y,z 

Partial derj 

M.C 

Partial deri 

00 

Infinity 


respect to x,y,z 
respect to 5,n,t 
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COMPUTATIONAL METHOD 


Flow Equati on 

The classical transonic small -disturbance flow equation written in terms of 
the disturbance velocity potential has been extended by retaining additional 
terms found in the full potential equation. These terms make it possible to 
resolve shock waves with appreciable sweep on wings in the X-Y (wing) plane. 
The global flow equation is: 


+ [l - (y 



( 1 ) 


This equation is also used in each component embedded grid system. The only 
exception is the winglet embedded fine mesh system. Wing-tip-mounted verti- 
cal wmglets may, at certain conditions, exhibit shock waves, swept in the X-Z 
(winglet) plane. To improve the ability to capture this type of shock wave, 
the following flow equation is used: 


1 - M 


2 


•xi 


(y + Dm 2 ® 

oj x 


Y + 1 
2 



2 

2 Mo 

.O Z 


0 

xz 


(2) 

+ [l - (y - 1 )M 2 ® j, +3 =0 

x zz yy 

In Equation (2), the X-Y cross-flow terms have been replaced by equivalent 
X-Z cross -flow terms. 


Pylon- type surfaces, in general, are not highly loaded. Shock waves are 
expected to propagate normally off neighboring components. Since swept pylon 
shock waves are not expected in most applications, special provisions involving 
the flow equation have not been made. 
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Computational Grid Approach 

The feature which distinguishes this approach from that of other methods 
is the use of multiple embedded grid systems. Conventional schemes employ a 
single continuous computing grid system. The multiple grid approach provides 
many advantages, if the analysis of complex shapes is the primary concern. 
Figure 2 illustrates a transport configuration with the aircraft component grid 
system boundaries included. The entire arrangement is positioned in the 
center of a global Cartesian coordinate system. The global grid is stretched 
so that boundaries represent infinity. This grid also serves a second pur- 
pose; it provides a means for interactions between the component embedded 
grid systems. Crude and fine grid systems interact or communicate by using 
an overlap region wherein the flow is computed twice for each iteration. This 
procedure is described in Reference 6. 



R84-1 137-002 D 


Figure 2 Multiple Grid Approach for Complex Aircraft Flow Simulations 
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The embedded grid systems provide computational resolution that is much 
more detailed than that of the surrounding global crude mesh. A detailed 
analysis is performed only in a region close to the configuration where gradi- 
ents are high. Ideally, three-dimensional methods should permit grid densities 
that are higher than two-dimensional and axisymmetric counterparts. This 
need is generated by the increased complexity of a flow field which has an ad- 
ditional dimension. The present method fulfills this requirement. There are 
100 evenly spaced mesh points available in the streamwise direction along fuse- 
lage, nacelle and pod surfaces. In addition, wing, pylon and winglet chords 
are modeled using 100 streamwise points. 

Discretization errors at the crude/fine grid interface are a function of 
both mesh spacing and flow gradient. Fine grid boundaries are positioned 
sufficiently far from the boundary surface to insure that these errors are 
small. Numerical experiments were used to "set" grid interface boundaries. 

Pod and Engine Nacelle Surfa ces 

The conventional technique generally applied to body or fuselage shapes 
was initially used for simulating pod and nacelle surfaces. This required that 
a computational surface representing the nacelle extend from upstream to down- 
stream infinity. Nacelle boundary conditions were applied along the length of 
tube starting at the inlet face and ending at the exhaust exit. Computations 
indicated that this approach would not be satisfactory. It is suspected that 
the problem is related to the difference between the physical and computational 
streamtube surfaces (Figure 3). The difficulty is compounded since these 
streamtubes typically pass very close to wing or body surfaces. Physical 
streamtubes are deflected by wing upwash/downwash fields. In addition, they 
are contoured by inlet mass flow ratio and exhaust pressure ratio effects. The 
computational streamtube is constrained to lie along existing grid lines. These 
tubes may be contoured by appropriate slope- type boundary conditions, but 
the physical position at which the boundary condition is applied is still con- 
strained to lie on the computational grid surface. This erroneous streamtube 
interference problem increases as wing-nacelle separation decreases. As might 
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be expected, the problem is more severe when streamtubes pass over a lifting 
wing at transonic conditions than it is when streamtubes pass below the wing. 

Exact position of inlet and exhaust streamtubes is difficult to determine. 
If position was known, or if it could be calculated with confidence, it could 
not be modeled easily. The best solution to this problem might be to eliminate 
the streamtube completely. Appropriate flow field potentials can be assigned 
to inlet and exhaust surfaces based on the inlet mass flow ratio and nozzle 
pressure ratio. The flow will then develop naturally in front of and behind the 
nacelle. This scheme should then also be suitable for modeling a pod or closed 
store which does not generate inlet or exhaust streamtubes. 



COMPUTATIONAL STREAMTUBE 



Figure 4 illustrates crude and embedded fine grid arrangement for a typi- 
cal nacelle or pod surface. Two operational modes are available: if nacelle or 

pod interference or nearby wing and fuselage surfaces is the primary concern, 
the crude grid representation is sufficient; if details of pressures on the 
nacelle surface are required, the fine grid solution can be obtained. 
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GRID SYSTEM GRID OVERLAP REGION 

R84-1 1 37 -004D 

Figure 4 Crude/Fine Grid Arrangement for Pods & Nacelles 

Computations indicate that a four-point diamond pattern provides the best 
approximation of the pod surface in the crude grid system. Grid points above 
and below and to the left and right of the pod center point make up this pat- 
tern. Note that this "best" flow simulation pattern may not be the same for 
other computational methods, because it depends on both the Y and Z mesh 
spacing used and the typical size of aircraft nacelles. Approximately 10 to 20 
mesh points fall along the nacelle in the axial direction. The actual number 
depends on the length of the nacelle relative to other aircraft dimensions. 

Approximately 28 mesh points represent the nacelle at each axial station if 
the embedded grid system is used. One hundred evenly spaced points fall be 
tween inlet and exhaust surfaces. The embedded fine grid interacts with the 
global crude grid in the same manner as the wing and fuselage grid systems. 
This requires a grid overlap region (Figure 4) . Crude grid flow field poten- 
tials are interpolated (linearly) to determine fine grid boundary potentials 
around the fine grid perimeter. Note that this perimeter also represents the 
outer boundary of the overlap region. The fine grid system is then relaxed to 
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determine the detailed flow between the outer boundary and nacelle surface. 
Resulting fine grid field potentials are interpolated to update crude grid points 
which make up the inner boundary of the overlap region. The global crude 
grid system is then relaxed and the cycle is repeated. The solution process 
continues by alternately sweeping the fine and crude grid systems until both 
are satisfactorily converged. Detailed pressure distributions are obtained by 
differencing fine grid nacelle boundary potentials. 

Figure 5 illustrates a typical nacelle boundary point along with field 
points required for determining boundary potential values. For a nacelle 
shape defined by: 


f(X,Y,Z) = 0 (3) 

the small disturbance boundary condition is: 


f + f t> + t : = 0 

x y y 7 . ; ■. 


(4) 





V 

9 


, N, a, j3, 9 , AY, AZ 


FIELD' 


[V , Ar, 9 N ) 

\ c LT J / 


r - GEOM SURFACE RAD 
9 

r = COMP SURFACE RAD 
C 


Figure 5 Nacelle & Pod Side Boundary Points 
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Three -point extrapolated differences are substituted for the velocities in 
Equation (4) . The boundary surface representing the nacelle is constrained to 
lie along Cartesian grid lines. As a result, there are regions along the nacelle 
where the computational nacelle surface and the geometric or physical nacelle 
surface do not coincide. Slender body theory is used to provide a boundary 
condition correction for both lifting and non-lifting displacement effects (Ref- 
erence 6, page 28). The thickness correction is proportional to the ratio of 
the geometric and computational nacelle radii. The correction due to flow inci- 
dence is proportional to the ratio of the two cross-sectional areas. The nacelle 
boundary potential value is: 


where 


and 


$ c ( j , k.) = A/B 


(5) 


A = N + N 

x r ' y 


4 $ 


Li * U) ' *(i ♦ 2.k) M 


2AY 


pod k r 


+ N 


4<t> 


(j> k * 1) Yi.k + 2) 

2AZ ( “ + “pod^T^ 

C 


(5.1) 


B =| 


N N 
+ -1 
AY AZ 


(5.2) 


When the solution process is complete, computational surface velocity dis- 
tributions must be converted to provide velocity distributions on the true 

nacelle surface. Once again, slender body theory provides the required cor- 
rection : 


v = V + (r - r ) Is (li) 

8 c g c dx r c ' 3n ' pod 


( 6 ) 
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Since inlet and exhaust streamtubes are not used in the present ap- 
proach, it is necessary to develop boundary conditions for inlet and exhaust 
surfaces. An effort has been made to minimize the number of variables re- 
quired for specification of fore and aft flow conditions. 


First, consider the exhaust surface. By using the energy equation, an 
expression can be written for exit velocity (Reference 7). 


V 


J 



(7) 


A similar expression can be written for the freestream velocity at infinity. 
The following expression evolves for the jet to freestream velocity ratio: 



Note that if system entropy losses are assumed to be small, jet 
static pressure are approximately equal; i.e.. 


( 8 ) 


and freestream 


P 


J 


s 


p 

“s 


(9) 


and nozzle pressure ratio is defined by: 


NPR = P , /P 

J 0 S (10) 


In addition, Reference 8 provides an expression for freestream static to total 
pressure ratio: 


P 

n 

s 


p 

.tl 

0 



Too 

v - 1 


(ID 
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For a freestream velocity that is equal to unity, the expression for jet exit 
velocity becomes : 


v 


J 





1 - (l/NPR) Y J 


y - 1 7 

1 - 1/(1 


( 12 ) 


For simulating exhaust effects typically found in experimental data, the jet ex- 
haust is cold and the following relations can be assumed: 


Y, = Y =1.4 
J ™ 


(13.1) 


Equation (12) reduces to: 



(13.2) 


1 - (l/NPR) 


0.2857 


cold 


l 


(l/ (1 + 0.2FT)] 


(14) 


For simulating flight test exhaust effects which typically involve hot jets, the 
following assumptions are made: 


(15.1) 

(15.2) 

(15.3) 


Equation (12) reduces to: 

V J h 

hot 


2.0636 


(l/NPR) 


0.2481 


T J = 1500°F + 460 J F 
T = 70°F + 460°F 

CD 

Y J = 1.33 


1 


2 


Figure 6 illustrates an engine nacelle with grid points that are involved in 
computing inlet and exhaust face boundary potentials. By using a three-point 
extrapolated difference operator, an expression for exhaust velocity can be 
written: 


V =l + o 
e x 


1 + 


1 

2 AX 


4<t> 


(I 


1) 


3 *(I ) " 5 (I 
e ( 


+ 2 ) 


(17) 


Z FREESTREAM VELOCITY, 
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Figure 6 Nacelle Inlet & Exhaust Boundary Points 

Exit velocity is then set to be equal to jet velocity providing the following 
expression for exhaust flow field potential: 


t> 


I 

e 


2(1 - V )AX + 4d) (I + 

e 

3 


1) $ (I +2) 

e 


(18) 


Values of V, will depend on whether a cold or hot jet is being simulated. 

A similar expression can be derived for the value of the flow field poten- 
tial at the inlet face. In this case, inlet mass flow ratio provides the required 
velocity relation: 
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MFR 


(19.1) 


V 
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r 


v_ = MFR = 1 + d) 

1 x 


(19.2) 


or 


MFR 


1 

2 AX 




(i. ) 

i 


— U ■ 4 ' 


(1. - 
1 


1) 


' (I ! - 2) 


The inlet flow field potential becomes: 

2 ( MFR - 1 ) AX + 4$ 


1) ®(I. - 2) 

I 


( 20 ) 


( 21 ) 


For pods and stores, or other bodies which do not have inlet and exhaust flow 
fields associated with them, both inlet and exhaust velocities are set to 0 
(MFR = 0), to simulate stagnation conditions. Tail and nose flow field potenti- 
als are then computed according to Equations (18) and (21), respectively. 


Positioning of the computational nacelle surface in the global crude grid 
structure is an important aspect of obtaining a good flow simulation. The ver- 
tical or Z grid system is constructed using a tangent function for stretching 
and the wing average chord length as a characteristic length (see Reference 6, 
page 12). If existing mesh points are interrogated to find the closest point to 
the nacelle centerline, more often than not, a poor flow simulation results. 
Calculations indicate that computed interference effects between the nacelle and 
wing are very sensitive to small changes in nacelle position. For this reason, 
it becomes very important to assure that computational and physical nacelle’ 
heights are identical. For configurations with nacelles or pods, the basic ver- 
tical grid system is adjusted (compressed) until the grid point just beyond the 
nacelle centerline position (as measured from the wing plane) falls on that po- 
sition. This grid adjustment procedure has been illustrated in Figure 7. 
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Figure 7 Global Crude Grid Adjustment for Nacelle Positioning 


Span wise positioning 1 of the nacelle surface is not as critical as vertical 
positioning in most applications. The existing spanwise grid system places 18 
mesh lines between the centerline and wing tip. Positioning the wing tip be- 
tween two mesh lines to simulate proper aspect ratio is still the most critical 
requirement. As a result, the computational nacelle center may be positioned 
as much as 3% half- span away from the physical nacelle center position. This 
type of discrepancy will, as a rule, not lead to simulation problems. But in 
certain applications (i.e., fuselage mounted nacelle), the spanwise nacelle posi- 
tion may be as critical as the vertical position. Special grid provisions or 
compromises might be required to successfully simulate this type of flow. 


Pylon Surfaces 

Wing surfaces are modeled in an embedded grid system that is skewed and 
tapered to "fit" the wing planform (Figure 2). The spanwise lines (n) of this 
grid system are aligned with wing constant percent chord lines. As a result, 
a smooth simulation of wing leading edge flow is obtained and chordwise reso- 
lution is uniform between the root and tip of the wing. 
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Modeling of a surface which is not aligned with this grid system presents 
a problem. A symmetry plane and fuselage side are good examples of nonalign- 
ed surfaces (see Reference 6, page 24). Numerical stability is impaired in a 
region near the nonaligned surfaces and, in some cases, the relaxation solution 
diverges. Some investigators have unswept grid lines in these critical 
regions; however, the resulting local Cartesian grid typically produces erro- 
neous oscillating pressure distributions. Pylon surfaces, like fuselage sides 
and symmetry planes, are not aligned with existing wing grid structure. 
Thus, the primary problem that must be overcome if pylon surfaces are to be 

modeled, is due to the inherent numerical instability of the pylon boundary 
condition . 

Existing grid lines which cross the wing planform in the spanwise direc- 
tion are generated by a hyperbolic tangent function and the requirement that 
the wing tip be centered between two grid lines. Eighteen lines are con- 
strained to lie between the symmetry plane and wing tip. The computational 
pylon surface is positioned at one of these existing spanwise grid lines. A 
discrepancy which may be as large as 3% of wing semi-span can exist between 
the physical and computational pylon span position. Note that this is also true 

for nacelle positioning. No attempt has been made to correct for this displace- 
ment effect. 


Figure 8 illustrates arrays of mesh points which are embedded along the 
wing. At any particular span station, the X-Z system is Cartesian in charac- 
ter. This wing grid system is also used to represent the pylon boundary sur- 
face. A separate grid system specifically for the pylon surface is not used. 
In addition, special provisions for swept leading edge simulation in the rectan- 
gular grid have not been made. As a result, current pylon modeling permits 
camber, twist and thickness variations; but the pylon planform is constrained 
to have no sweep or taper. This limitation has been sketched in Figure 9. 
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Figure 8 Pylon Surface Boundary Points in Wing Fine Grid System 


PHYSICAL PYLON SURFACE 
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Figure 9 Physical 8t Computational Pylon Surfaces 






The standard wing fine grid system vertical boundaries are positioned 
based on a percentage of wing root chord length. When pylons are modeled, 
however, the fine wing grid lower boundary must be adjusted to insure that 
the pylon lower surface does not project below the grid limit. This is accom- 
plished by positioning the lower grid boundary at the pylon base level and al- 
tering the vertical grid spacing so that the total number of vertical grid points 
remains within original limits. 


There are several approaches that can be used to enforce pylon boundary 
conditions. In this effort, an attempt has been made to find the best compro- 
mise between providing a good flow simulation and providing a stable numerical 
solution that will converge for a variety of different flow/geometry cases. 
Calculations made to date indicate that if extrapolated differences are used to 
’set” the inboard and outboard pylon surface potential values, instabilities in 
the pylon leading edge region will result. Similarly, if pylon boundary con- 
ditions are used as a constraint on the cross-flow velocity, the solution will 
diverge. By using dummy potential arrays, however, a scheme can be devel- 
oped which actually enhances the numerical stability of the solution. Calcula- 
tions indicate that a satisfactory flow simulation is also obtained. 

Figure 10 illustrates which grid points are involved in generating dummy 
potential values for inboard and outboard pylon surfaces. In the skewed grid 
system, the velocity component, v, can be written: 

v = 4> = $ n + 

y n y i y 

This can be set equal to the pylon boundary slope: 

P y f pyi S pyi 

wkere the sign convection is due to the inboard surfaces of both winglets and 
pylons being specified as the input airfoil section "upper” surface. A 
special two-point (first order accurate) difference operator can be used to 
provide an approximation of the ^ and ^ derivatives. The expression for 
the inboard surface dunmy potential point becomes: 
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A similar expression can be written for the outboard surface dummy 
point: 
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For each side, note that the coefficient of the dummy potential point is greater 
than the coefficient of other points involved. This enhances the effective dia- 
gonal dominance of the system. 



R84-1 137-010D 

Figure 10 Grid Points Required for Computing Pylon 
Inboard/Outboard Boundary Conditions 
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The computational method operates with three planes of grid point poten- 
tials in core at one time. When a pylon station is reached, the outboard plane 
is modified in the region which is directly aligned with the pylon surface. 
Equation 24 provides the temporary dummy flow field potential values. The 
inboard pylon surface is then relaxed in a manner similar to that used for 
normal field points. When the wing/pylon trailing edge is reached, the origi- 
nal outboard potential values are replaced, and inboard plane dummy values 
are established. The outboard side of the pylon surface is then relaxed and 
the pylon surface simulation is complete. Vertical pylon circulation is com- 
puted by: 


r(z) 


pyl 
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This value of circulation provides the jump condition required for differencing 
across the pylon wake surface: 
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Winglet Surfaces 

For wing- tip-mounted winglets, the small size of the lifting surface alone 
presents a formidable obstacle. Winglet planform areas are typically between 
1/40 and 1/70 the area of the main wing. Sufficient resolution must be pro- 
vided for both wing and winglet surfaces simultaneously. The stability prob- 
lems which occurred when pylon development was in progress are not incurred 
since rectangular grid systems are set up specifically for winglet surface sim- 
ulation to overcome the resolution problem. 
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Spanwise grid arrangement for a conventional wing tip is shown in Figure 
11. The wing tip is positioned midway between two grid lines. This provides 
the best approximation of wing tip aerodynamics. If a winglet is modeled, the 
grid system must be adjusted to provide grid lines at the wing tip. These 
wing tip points are used for both the winglet boundary surface and the ver- 
tical winglet wake which extends to downstream infinity. 
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Figure 11 Wing Tip Spanwise Grid System for Winglet 
& Conventional Type Planforms 


If a fine grid system were embedded along the winglet surface in the 
same manner currently used for wing surfaces, a diverging solution would 
probably result. This would be caused by the large disparity in resolution 
between the existing global crude grid structure and the embedded fine grid 
system. To overcome this problem, an intermediate mesh system is placed 
around the wing tip region. Figure 12 illustrates the resulting three-level 
grid arrangement required for resolving details of winglet flow fields. Both 
the global crude grid and intermediate grid are rectangular in character. The 
winglet embedded (fine) grid system is skewed and tapered to fit the winglet 
planform. 

During the first phase of the solution process, only crude and intermed- 
iate grid systems are active. Winglet boundary conditions are imposed in the 
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higher resolution intermediate grid system. Since vertical line relaxation is 
used, the boundary condition is enforced by specifying the cross-flow velocity 
component in Equation (1). (Note that this will involve both the 6 and <b 
terms.) Winglet circulation and wake jump computations are identical to those 
specified for pylon surfaces (Equations 26, 27 and 28). Crude grid winglet 
points form a fixed potential surface and flow field potentials are generated by 
the interpolating values found in the intermediate grid. Most important, how- 
ever, the crude grid system provides a means for carrying vertical wake cir- 
culation downstream. 



GLOBAL 
CRUDE GRID 


Figure 12 Crude/Medium/Fine Grid Arrangement for 
Wing-Tip-Mounted Winglets 


The winglet flow field computed during the first phase of the solution 
process is used to provide a starting point for the second phase which re- 
quires crude-medium-fine grid interactions. Fine winglet grid boundary con 
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ditions and fine wing grid boundary conditions are identical. Only the coordi- 
nate direction is changed. For the winglet inboard surface, the boundary con- 
dition is: 
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The outboard surface condition is: 
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During the second phase of the solution process, flow field potentials are fixed 
on both crude and intermediate grid winglet surfaces. The global crude grid 
solution provides outer boundary potentials for the intermediate grid. The in- 
termediate or medium grid solution provides outer boundary potentials for the 
embedded fine grid system. 


Solution Process 

This section describes the procedure which is used to solve for the multi- 
ple wing and body surface flow field. The procedure includes steps that are 
executed when nacelles, pylons and winglets are attached to the basic wing- 
body configuration. For simpler configurations, individual component steps are 
simply bypassed. Figures 2, 4 and 12 may prove to be useful in visualizing 
operations which are to be described. The solution process can be broken 
down into two separate phases. 

Phase 1: INITIAL CRUDE GRID FLOW FIELD SOLUTION 

An initial solution is obtained with the wing-body represented in 
the crude mesh. Typically, 100 iterations are sufficient. After 60°o 
of the crude cycles are completed, flow field potentials from the 
global crude grid are used to establish a starting flow in the wing 
tip augmented grid (TAG). Note that TAG was referred to as an 
intermediate grid in the winglet surface section. For the remaining 



40°6 of the crude grid cycles, winglet boundary conditions are en- 
forced in TAG and the global grid serves also to carry wing and 
winglet wake systems downstream. 


Phase 2: 
Step 1: 


Step 2: 

Step 3: 
Step 4: 
Step 5: 
Step 6: 


CRUDE/FINE GRID INTERACTION 

Using global crude grid flow field potentials obtained in Phase 1, 
a starting flow field for the embedded wing grid is established. A 
similar initialization procedure is performed for the embedded pod or 
body grid system. Flow field potentials from TAG provide a start- 
ing flow field for the winglet embedded grid system. Pylon out- 
board potential arrays are set equal to inboard potential values 
which make up the basic wing flow field. Pylon circulation distribu- 
tions are set to zero. 

The wing fine grid system is relaxed with outer perimeter poten- 
tial values fixed and conventional Neumann boundary conditions at 
fine grid section boundary points. Crude grid wing boundary point 
potentials are updated based on results from this step. 

The winglet fine grid system is relaxed in a manner similar to 
that in Step 2. Both the crude grid and TAG winglet surface 
boundary potentials are updated. 

Wing boundary layer 6* is computed. Embedded body fine grid 
perimeter points are updated, based on global crude grid field po- 
tentials. Body fine grid boundary potentials are computed. 

The body fine grid system is relaxed subject to fixed potential 
values on the grid perimeter, body surface, and wing plane surface 
(see Figure 15, Reference 6). 

The crude grid body surface is updated based on fine grid field 
potentials determined in Step 5. 
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Step 7 : 


The nacelle fine grid system is relaxed in a manner similar to 
that of the body (Steps 5 and 6). In addition, however, flow field 
potentials at inlet and exhaust faces must be determined before each 
relaxation sweep. 

Step 8: The global crude grid is relaxed using Dirichlet boundary poten- 

tials generated by Steps 1-7. Field potentials from the crude grid 
now provide outer boundary potentials for embedded grid arrays. 

Step 9: The TAG is relaxed with fixed perimeter and winglet surface 

boundary points. 

Steps 1 through 9 are repeated until the grid systems are satisfactorily 
converged. Typically, 80 cycles are sufficient. When the boundary layer op- 
tion is selected, a boundary layer displacement thickness (6*) will be com- 
puted every 20th iteration starting with the fifth iteration in phase two. The 
6* gradient is added to the wing surface slope for the wing boundary condi- 
tion required in step 2. If the fine body or pod grid option is selected, grid 
perimeter potentials arc updated every 10th iteration and the grid systems are 
relaxed every second iteration to conserve computing resources. 

Pressure, F orce and Moment Coefficien ts 

Pressure, force and moment coefficient calculations for wing and body 
components are described in Reference 6. All pressure coefficients on pod, 
pylon, and wingiet components are computed in the same manner as those de- 
scribed in Reference 6. Pylon surfaces, which are perpendicular to the wing 
plane, do not contribute to the total configuration forces and moments. For 
this reason, only nacelle and winglet force and moment calculations will be de- 
scribed . 



Winglet section coefficients are obtained by integrating pressure coeffi- 
cients on both inboard and outboard surfaces. 
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In addition, a sectional moment about the axial reference position is computed. 
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Winglet lift, moment and drag coefficients are obtained by integrating sec- 
tion coefficients across the winglet span and multiplying by winglet cant angle. 
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Winglet skin friction drag is computed by using the Prandtl-Schlichting formula 
corrected for compressibility effects (see Equation 64, Reference 6). 

Pod or nacelle force and moment coefficients are obtained by integrating 
surface pressure distributions and adding a skin friction component based on 
wetted area. A pod cross-sectional force coefficient distribution is computed 
first. 
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The longitudinal coefficients are then integrated along the length of the body. 
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As was the case for winglets, the pod skin friction coefficient is computed 
using the Prandtl-Schlichting formula (see Equation 64, Reference 6). 



The total configuration force and moment coefficients are obtained by add- 
ing the various component coefficients. 
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Two-Dimensional Airfoil Analy sis 

The modified small disturbance equation/planar boundary condition ap- 
proach, the basis for the present method, may in some cases produce results 
which are inferior to those that would be obtained if a full potential 
equation/exact boundary condition method were available. These discrepancies 
might become significant when surface shapes have extraordinarily high gradi- 
ents or when flow conditions are extreme. Unfortunately, it would be difficult 
to evaluate potential problem areas when complex configurations are of 
interest. Flow interference can be complex and many different components of 
the computational method are contributing to the final result. For these 
reasons, special provisions have been made to permit analysis of two- 
dimensional airfoil shapes. Abundant experimental data and easy-to-use full 
potential equation/exact boundary condition computer analyses can then be 
used to assess problem areas in two dimensions before complex three- 
dimensional analyses are performed. 

Figure 2 illustrated the position of the wing embedded fine grid system. 
This mesh array is made up of a series of two-dimensional arrays (see Figure 
8) . If a single planar array is used at the symmetry plane, a two-dimensional 
airfoil can be treated. Unlike the wing-body case (see Figure 13), however, 
far field boundaries which represent infinity are non- zero since the airfoil ap- 
pears as a concentrated point vortex. This has been illustrated in Figure 14. 
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Figure 13 Boundary Conditions at Infinity for Wing-Body Analysis 
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Figure 14 Boundary Conditions at Infinity for Airfoil Analysis 



The global flow equation (Equation 1) which is used for all three- 
dimensional computations is simplified by eliminating all cross -flow terms when 
airfoil analysis is required. 
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Airfoil pressure coefficients are computed using the following equation. 
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COMPARISONS & TYPICAL RESULTS 

Correlations with experimental data for isolated bodies, isolated wing’s, 
and wing- body combinations were presented in Reference 6 using the basic 
wing-body code. These examples will not be included in this report, but it 
should be noted that the present method will reproduce those results in ad- 
dition to providing the more complex configuration flow simulations contained 
herein. Thus, the present method can be used to analyze 1) two-dimensional 
airfoils, 2) isolated wings, 3) isolated bodies, 4) wing-body combinations, and 
5) aircraft configurations with as many as four pods/nacelles, four pylons, and 
wing-tip-mounted winglets. This provides flexibility for studying component 
interference effects. 

The computations which follow illustrate the new airfoil and nacelle, 
pylon, winglet simulation capabilities. Component interference effects are 
shown for cases where component on/off wind tunnel data is available. For 
some components (nacelles and winglets), detailed pressure comparisons are 
included. For one case (G-III), flight data is used to augment wind tunnel 
results . 

Airfoil Correlations 

Two airfoils have been selected for correlation studies. The NACA 0012 
airfoil is somewhat conventional in shape. It is also a section which is often 
used in theoretical or analytical test samples. The second airfoil selected is 
the new type NASA LS(1) airfoil which features considerable aft-loading. Ex- 
perimental data for these sections can be found in Reference 9 along with 
Korn-Garabedian and Carlson 2-D code analysis results. Correlations with ex- 
perimental data using the present method can be seen in Figure 15. 
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Figure 15 Two-Dimensional Airfoil Pressure Distribution Correlation 



Nacelle Flow Simulations 


Two different transport configurations have been selected for comparisons 
to experimental data. Together, they illustrate both nacelle operational modes. 
The first, a crude grid nacelle simulation, provides flow disturbances on the 
wing. The second involves a fine embedded grid computation which is used to 
predict detailed pressures on the engine nacelle surface. 

Figure 16 illustrates the G-III configuration. This transport has a low 
wing and two fuselage-mounted nacelles. The wing has an aspect ratio of 6, 
and a taper ratio of 0.279. Wing leading edge sweep is 31.7 degrees. 
Nacelles are canted 2.75 degrees nose up, and 0.5 degrees nose inboard. 




Figure 16 G-III with Pressure Tap Locations for Nacelle Interference Study 

G-III cruise conditions occur near Mach 0.75 and an angle- of -attack of 4 
degrees. The wing root pressure distribution at this flow condition is rela- 
tively shock free and not very challenging for the transonic computational 
method. For this reason, a higher Mach number case has been selected for 
comparisons. The Mach 0.85 flow condition exhibits a strong shock wave. 
The selected incidence must be lower than cruise levels to avoid flow separa- 
tion. 
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Wing upper surface, superimposed, computed pressure distributions for 
the basic wing-fuselage and the wing-fuselage-nacelle cases can be seen in 
Figure 17 at Mach 0.85 and 0 degree angle-of-attack. A single shock wave 
spans the wing between the fuselage side and wing tip. With the nacelle 
mounted, the shock wave disappears in the wing fuselage- nacelle juncture re- 
gion. Resulting flow on the wing under the nacelle exhibits a deceleration in 
front of the inlet face, and an acceleration just behind the inlet lip. The inlet 
MFR for this case is 0.66. 



Correlations with experimental data at three span stations along the wing 
can be seen in Figure 18. Note that the character of this complex 
three-dimensional flow field is predicted quite well. The lower surface pres- 
sures are essentially unaffected by the presence of the engine, and the 
interference effect can be observed to decay properly as distance from the 
nacelle increases. Experimental data is taken from Reference 10. 
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Figure 18 G-III Nacelle On/Off Wing Pressure Distribution Correlations 
M = 0.85, a = 0° 



The C-5A configuration is illustrated in Figure 19. For flow conditions 
near the cruise point (M^ 0.775 and a = 2.0 degrees), an embedded fine grid 
has been positioned around the inboard engine nacelle. Experimental data can 
be found in Reference 11. The actual engine consists of a fan cowl, turbine 
cowl and plug. Engine surface discontinuities and expected flow separation aft 
of the fan cowl cannot be properly modeled using the present attached flow 
methodology. For this reason, the long duct flow- through nacelle (Configura- 
tion 10, Reference 11) is modeled for analysis. Computed results are compared 
to experimental pressures which were measured on the basic nacelle (Config 
uration 1 ) . 




Correlations with experimental data at three stations around the nacelle 
can be seen in Figure 20. Note that both basic and flow-through nacelle 

shapes have been sketched. Computed results on the turbine portion of the 
nacelle can only be considered as approximate. They are included because the 
character of this flow which is influenced by the wing leading edge is properly 
predicted. 

Figure 21 shows inboard nacelle port and starboard comparisons super- 
imposed. Hie variation in flow character is probably caused by a combination 
of interference effects that include fuselage effects inboard, second nacelle 
effects outboard, and the swept back wing leading edge. 
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Figure 21 Comparison of C-5A Inboard Nacelle 

Port and Starboard Pressure Distribution 
Correlations 

M = 0.775 a = 2° NPR = 2.84 


Pylon Interference Effects 

The C-5A configuration (Figure 19) has been selected to demonstrate the 
pylon computational capability. Reference 11 provides experimental data for 
the basic C-5A wing-fuselage configuration and various configurations with dif- 
ferent pylon and nacelle combinations. Wing pressure data which illustrate 
pylon interference effects are available at both pylon span stations (Figure 
22). For the wing lower surface, pressure taps are positioned on both sides 
of the wing pylon juncture. A double shock wave system exists inboard of the 
mid-semi-span region; outboard, only a single shock wave exists. Correlations 
with experimental data at the inboard station for the pylon/nacelle on and off 
can be seen in Figure 23. Outboard station comparisons can be found in 
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• WING UPPER/LOWER SURFACE 
PRESSURE TAP LOCATIONS 
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Figure 22 C-5A Wing Stations for Pylon/Nacelle Interference Study 
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Figure 23 C*5A Inboard Engine Station Pressure Distribution Correlation 
0.775 a = 2° NPR = 2.84 
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Figure 24. Note that pylon/nacelle interference pressures are properly pre- 
dicted and the multiple shock wave character is accurately reproduced. These 
comparisons were made for Configuration 1 in Reference 11. The P./P^ value 
of 1.5 converts to an NPR of 2.84. 
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Figure 24 C-5A Outboard Engine Station Pressure Distribution Correlation 
M = 0.775 a - 2° NPR = 2.84 

oo 


Winglet Flow Simula tions 

Reference 10 provides experimental data suitable for verifying the 
winglet computational capability. The G-III configuration with labeled experi- 
mental pressure tap stations can be seen in Figure 25. Figure 26 illustrates 
the wing tip comparison for both winglet on and off. It is interesting to note 
that the more complex winglet/ wing tip arrangement shows better agreement. 
This is probably caused by differences in wing tip rake for the two planforms. 
The analysis method cannot resolve details of the highly swept tip-leading-edge 
contour, thus, it predicts higher loadings in the leading edge region than 
those given by experiment. The winglet-on wing tip is trapezoidal in charac- 
ter and, therefore, in better agreement with the wing tip region computational 
model. 
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Figure 25 G-lll Wing/Winglet Pressure Tap Locations 
for Evaluating Winglet Interference 
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WINGLETS ON 

SECTION CL = 0.505, CM = -0.01 7, CD = 0.0057 


Figure 26 G-lll Wing Tip Pressure Distribution Correlation 
M = 0.75 a = 4° 




Winglet pressure distribution correlations at three stations along the 
winglet span can be seen in Figure 27. Viscous effects are not modeled in 
these computations. A mild shock wave can be identified as it follows the 
leading edge of the winglet on the inboard side. 
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Figure 27 G-lll Winglet Pressure Distribution Correlations 0.75 



SECTION CL - 0.580, 
CM - -0.091, CD = -0.01 14 


a =4° 


Combined Interf erence Effects 

Pressure distributions across the G-III span show the good results that 
can be obtained after modeling all configuration components (see Figure 28). 
This comparison includes both wind tunnel and flight data. 

The KC-135 provides an additional case that can be used for code valida- 
tion (Reference 12). This configuration has four pylons, four engine nacelles, 
and wing- tip-mounted winglets. Pressure data is available at several wing 
stations and several winglet stations. Strong shock waves are not present at 
cruise conditions. If flow conditions are pushed beyond the cruise point, wing 
tip flow separation is apparent. Comparisons with experimental data at the 


48 


I 




Figure 28 G-lll Wing Pressure Distribution Correlation M = 0.78 a = 


cruise point can be seen in Figure 29. It should be pointed out that the 
winglet is positioned at a 90 degree angle to the wing surface when modeled in 
the computational method. In reality, both the KC-135 and G-III winglets are 
canted outboard 15 degrees. 

The C-141 transport configuration can be seen in Figure 30. This case is 
interesting because it illustrates interference effects caused by nacelles, 
pylons, and fuselage fairings. The fuselage geometry model, made up of body 
and cross-section lines, can be seen Figure 31. Wing-body juncture and land- 
ing gear fairings can be identified. 

Figure 32 compares unpublished wind tunnel pressure data taken at 
four wing stations. The flow condition is 0.77, a = 1.5 degrees, and R g 

2 X 10^. Note the analysis method angle-of-attack is 1/2 degree higher than 
that of the experiment. The experimental variations in flow field character 
which result from the addition of nacelle/pylon combinations are predicted quite 

well. 


*Wind t unn el data obtained from M. Imres, lockheed-Georgia Company. 
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CONCLUDING REMARKS 


The computational method described has been specifically developed to 
provide an engineering analysis for realistic aircraft configurations at transonic 
speeds. It is unique in its ability to treat a variety of shapes. This feature 
should prove to be useful in the study of aerodynamic interference effects. 
Compared to existing methodology, the approach provides very high computa- 
tional resolution. Resolution varies between 200 boundary points for a simple 

airfoil to over 9000 points for a complex wing-fuselage combination with 
nacelles, pylons and winglets. 

Method flexibility, which is required for treating a variety of complex 
shapes, is balanced by the simplicity of its components. Coordinate systems 
are essentially rectangular in character, a simple two-dimensional strip bound- 
ary layer analysis provides viscous corrections, and finally, a fast, easy to 
use fuselage modeling system yields arbitrary body shape surface normals. 
While more sophisticated components could be used, the simplicity and cost 
effective character of the present arrangement should enhance the probability 
of obtaining accurate flow simulations. 


55 


56 


I 



APPENDIX A 


COMPUTER CODE DESCRIPTION 


CONTENTS 

Page 

GENERAL COMPUTER CODE DESCRIPTION 59 

INPUT DATA FORMAT 64 

SAMPLE INPUT DATA SETS 73 

OUTPUT DATA (PRINT AND PLOT) FORMAT 86 

INPUT GEOMETRY VERIFICATION 92 

SUBROUTINE CALL SEQUENCE 97 

SUBROUTINE DESCRIPTION 102 

KEY VARIABLE DESCRIPTIONS HO 


^ecewto page punk not fmjwd 


57 


58 


I 



GENERAL COMPUTER CODE DESCRIPTION 

The computer code is operational on both IBM and CDC type com puters t 
Overlay structures are not used although this approach (for reducing' core re- 
quirements) may be advantageous depending on facility charging algorithms. 
The IBM version using the extended H compiler (opt = 2) requires approxi- 
mately 970 Kjq for storage and execution. There is considerable use of tempo- 
rary disk storage units. Since interpolation and searching is required, a re- 
sult of the mesh embedding approach, it is useful to have planar potential ( 0 ) 
arrays separate and addressable. As a result, 15 different units are currently 
employed. The disk unit number and a description of contents are listed be- 
low. 

DE SC RIPTION 

Input data transferred to Unit 1 (formatted data) 
Quick-geometry problem diagnosis printed- output 
Global crude grid potential array 
Fine wing grid potential array 
Fine body/nacelle grid potential array 
Fine winglet grid potential array 

Crude grid wing upper/lower surface boundary conditions 
Fine wing grid upper/lower surface boundary conditions 
Fine wing grid x- coordinate array 
Fine wing grid section surface ordinates 
Crude grid body surface normal (direction cosines) 

Fine grid body surface normal (direction cosines) 

Fine wing grid shearing angles 
Wing and body pressure coefficient arrays 
Boundary layer displacement thickness slope 

* — : ; 

Ibis appendix describes the IBM code version whereas Appendix B gives only 
the modifications for the CDC code version. 

pp pmy Nfi PAGE PLANK NOT FH-MSB 


DISK NUMBER 
1 
8 

12 

13 

14 

15 
80 
81 
82 

83 

84 

85 

86 

87 

88 
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Computer running time will of course vary depending on the facility and 
the mode of operation or operating system. The absolute levels specified may 
be out dated shortly after they are specified. IBM 370/3081 running times are 
specified below, however, since the relative increment for various options will 
remain essentially steady, these increments will be useful for estimating the 

TIME ( C PU Minute s) 

5 (50 crude/50 fine iterations) 

14 (100 crude/80 fine iterations) 

7 (150 crude/150 fine iterations) 

15 (100 crude/80 fine iterations) 

15 (100 crude/80 fine iterations) 

15 (100 crude/80 fine iterations) 

20 (100 crude/80 fine iterations) 

1 (No iterations) 

An effort has been made to minimize the amount of data required to define 
the configuration geometry and flow condition. This should simplify matters 
for most applications involving configuration analysis and reduce the chances 
for input errors. For example, the computational grid systems (extent and 
density) have been set in the FORTRAN coding to provide good results under 
most conditions. Occasionally, it will be advantageous to manipulate the preset 
values and limiters. FORTRAN coding changes will be necessary if this is the 
case. The following values and limiters may be modified in certain special ap- 
plications : 

1) Gas constant (5f = 1.4) 

2) Fine wing/body embedded grid limits or extent 

3) Fine wing/body embedded grid density 

4) Subsonic relaxation factor (w = 1.5) 

5) Boundary layer transition (X/C tran = 0.05) 

6) The number of in viscid cycles between each viscous calculation (cur- 
rently set to 20) 
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time and cost of using different options. 
CASE 

Isolated body 
Isolated wing 
Airfoil 

Isolated wing w/ viscous interaction 
Wing -body (body modeled in 
crude grid) 

Wing body w/pod- pylon -winglet 
Wing-body (body modeled in 
fine grid) 

Geometry/Grid verification 



All sample cases were computed using- the basic code without modification. If 
code modifications of type (2) or (3) are made care must be taken to insure 
that common and dimensioned arrays are sufficient and consistent. 

Input data format description can be found on the following pages. De- 
scriptions are thought to be relatively straightforward except in the case of 
wing section definition. It is important to extend wing planform/section defi- 
nition to the symmetry plane even for wing-body configurations. This serves 
several purposes. First, the code will compute a wing-body juncture which 
will be a function of both configuration geometry and the computational grid 
system. If the computational juncture is slightly inboard of the geometric 
juncture, section definition in this region becomes important. Second, the in- 
put planform shape provides both aspect ratio for the lift-induced drag com- 
putation and leference lengths and areas used to reduce integrated pressures 
to give force and moment coefficients. Finally, data input for defining a 
wing body configuration can be used directly for the isolated wing case. This 
feature can be used to study wing-body interference effects. 

It is important to note that provisions have been made for inputting a 
wing reference area. This number is used to reduce integrated pressure coef- 
ficients to provide wing or wing-fuselage total lift, drag, and moment coeffi- 
cients. All other reference areas and lengths are computed from input geom- 
etry and printed at the end of the output stream. If reference values used 
to reduce experimental data are different than those computed by the code, 
then computed force and moment coefficients must be rescaled. 

The computer code has been structured to permit the analyses of common 
aircraft configurations. There are a number of restrictions, however, that 
should be identified before modeling is started: 

1) No boundary layer 6* computed for winglet surfaces. 

2) No differencing approximations for secondary lower surface winglet. 

3) No provisions made for pylon surface attached to top of wing. 

4) The fine embedded nacelle/pod grid option cannot be used along with 
the pylon option; crude pod and fine pylon representation only. 
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5) The winglet cant angle is used only to resolve force and moment coef- 
ficients. Computational winglet surfaces are always perpendicular to 
the wing surface. 

6) The fine body grid and fine nacelle/pod grid potential arrays use the 
same dimensioned space . The two options cannot be implemented si 
multaneously . 

7) In addition to the current rectangular pylon planform restriction, the 
computational pylon lower surface is fixed at the wing fine grid lower 
boundary . 

8) Fuselage shoulder mounted nacelles must be positioned to allow at 
least one mesh cell between the nacelle and the fuselage. 


An overall view of the input data sequence can be seen in Figure A-l. 
Figure A-2 illustrates relationships between the wing to body distance (ZWING) 
and the wing to pod distance (ZBODY1). 
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Fig. A-1 Input Data Sequence 







INPUT DATA FORMAT 

jfc 

Excluding literal cards, all input data cards are punched in seven field 
ten digit format (7F10.0). A decimal point is required in each field. 


CARD 
NUMBER 
Card 1-A 


CARD VARIABLE 

COLUMN N AME DESCRIP TION 

1-80 TITLE Configuration or run title to identify graphic 

and printed output. 


Card 2-A 


1-10 


11-20 

21-30 

31-40 

41-50 

51-60 

61-70 


CASE CASE = 1. Isolated Body (omit cards -W) 

CASE = 2. Isolated Wing (omit cards -B) 
CASE = 3. Wing- Body 
AMACH Mach Number (AMACH<1.0) 

AOA Angle- of- Attack (degrees) 

g 

RE Reynolds Number (X10 ) 

AXIT Number of initial crude grid iterations. 

(AXIT = 0. for geometry verification only) 
AXITF Number of crude/fine grid iteration cycles. 

VISMOD VISMOD = 1. No viscous effects. 

VISMOD = 2. Viscous effects computed at 
end of in viscid analysis. 
VISMOD = 3. Inviscid/viscous interaction. 


NOTE: 

Omit card 3-A for CASE = 1 or CASE = 2. 

Card 3- A 

1-10 

PY 

Number of pylons(0 < PY < 2.) 


11-20 

VER 

Winglet code VER = 0. No winglet. 




VER = 1. Winglet definition to 




follow on cards WLT. 


21-30 

POD 

Number of pods or nacelles(0 < POD 5 2.) 


NOTE: Two-dimensional airfoils can be analyzed by setting ASECT - 1 
for the isolated wing CASE = 2. 


^Formatted input data files are created as card images. 
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CARD CARD 

VARIABLE 

NUMBER COLUMN NAME 

description 

NOTE: Omit card set W for CASE = 1 (card 2-A). 

Card 1-W 1-10 

ASECT 

Number of streamwise sections defining wing 
planform (1 < ASECT < 20). 

11-20 

ANIN 

Number of ordinates defining each wing sec- 
tion (ANIN < 60). 

21-30 

ANOSW 

ANOSW = 0. Sharp nose wing sections. 
ANOSW = 1. Blunt nose wing sections. 

31-40 

XMOM 

X-position about which moment is to be com- 
puted. 

41-50 

ZWING 

Z-position of wing (waterline), see Fig. A-2 

51-60 

REFAR 

Wing reference area. 

61-70 

WS 

Wing Cp distribution plot scaling per inch 
(typically 0.4 or 0.8). 

NOTE: Card set 

2-W through 5-W is repeated ASECT times. 

Card 2-W 1-10 

XPL 

Wing section leading edge (X-value). 

11-20 

YP 

Wing section span position (Y-value). First 
Y- value must be 0.0 (symmetry plane), even 
for wing-body case. 

21-30 

XPT 

Wing section trailing edge (X-value) . 

31-40 

TWIST 

Wing section local incidence (twist angle in 
degrees) . 

41-50 

AKODE 

AKODE = 0. Section ordinates identical to 
preceding section (omit cards 
3-W through 5-W) . 

AKODE = 1. New section definition expected 
on cards 4-W and 5-W. 

51-60 

AQXFG 

Number of fine X-grid points at wing tip 
station (default — 100). Read only on 
first Card 2-W. 
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CARD 

CARD 

VARIABLE 


NUMBER 

COLUMN 

NAME 

DESCRIPTION 

Card 3-W 

1-70 

XINW 

Wing section X coordinates (cards 3-W de- 
fined only for first wing section, ANIN 
values expected) . 

Card 4-W 

1-70 

YINU 

Wing section upper surface Y- coordinates 
(ANIN values). 

Card 5-W 

1-70 

YINL 

Wing section lower surface Y- coordinates 
(ANIN values). 

NOTE: 

Read Card 

set PY, PY 

times (Card 3-A). 

Card 1-PY 

1-80 

TITLPY 

Pylon title to identify printed output. 

Card 2-PY 

1-10 

PSEC 

Number of streamwise sections defining pylon 
planform (2 < PSEC < 10) . 


11-20 

PIN 

Number of ordinates defining each pylon sec- 
tion (PIN < 60). 


21-30 

PNOS 

PNOS = 0. Sharp nose pylon sections. 
PNOS = 1. Blunt nose pylon sections. 


31-40 

PSTA 

Y - position of pylon on wing. 

NOTE: 

Card set 

3-PY through 6-PY is repeated PSEC times. 

Card 3-PY 

1-10 

XPYL* 

Pylon section leading edge (X-value). 


11-20 

ZPYL* 

Pylon section height (Z-value). 


21-30 

XPYT* 

Pylon section trailing edge. 


31-40 

TPIST 

Pylon section local incidence in degrees (for 
positive TPIST - pylon nose rotation toward 
centerline) . 


*Note: Pylon planform description used for pylon wetted area calculation 

(pylon skin friction) only. Planform shape is now restricted to 
the boundaries of the wing fine grid system below the wing. 



CARD 

CARD 

VARIABLE 

NUMBER 

COLUMN 

NAME 

DESCRIPTION 

Card 3-PY 
(contd) 

41-50 

PKOD 

PKOD = 0. Section ordinates identical to 
preceding section (omit cards 
5-PY and 6-PY). 

PKOD = 1. New section definition to follow. 

Card 4-PY 

1-70 

XINP 

Pylon section X-ordinates (cards 4-PY de- 
fined only for first pylon section, PIN values 
expected) . 

Card 5-PY 

1-70 

YINUP 

Pylon section inboard surface coordinates 
(PIN values). 

Card 6-PY 

1-70 

YINLP 

Pylon section outboard surface coordinates 
(PIN values). 

NOTE: 

Omit card 

set WLT for WLT = 0 (Card 3-A). 

Card 1-WLT 

1-80 

TITLW 

Winglet title 

Card 2-WLT 

1-10 

VSEC 

Number of streamwise sections defining wing- 
let planform (2 < VSEC < 10). 


11-20 

VIN 

Number of ordinates defining each winglet 
section (VIN < 60). 


21-30 

VNOS 

VNOS = 0. Sharp nose winglet sections. 
VNOS = 1. Blunt nose winglet sections. 


31-40 

VANGL 

Winglet cant angle in degrees (for positive 
VANGL - winglet rotation upward from wing 
plane) . 
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CARD 

NUMBER 

NOTE: 

Card 3-WLT 


Card 4-WLT 

Card 5-WLT 
Card 6-WLT 

NOTE: 
Card 1-B 


CARD 

VARIABLE 


COLUMN 

NAME 

DESCRIPTION 

Card set 3 

i-WLT through 6-WLT is repeated VSEC times. 

1-10 

XVL 

Winglet section leading edge. 

11-20 

YV 

Winglet section span position (Note: Winglet 
is defined in plane of wing) . 

21-30 

XVT 

Winglet section trailing edge. 

31-40 

TVIST 

Winglet section local incidence in degrees 
(for positive TVIST, winglet nose rotated 
toward centerline) . 

41-50 

VKOD 

VKOD = 0. Section ordinates identical to 
preceding section (omit cards 
5-WLT to 6-WLT) . 

VKOD = 1. New section definition to follow. 

1-70 

XINV 

Winglet section X-ordinates (cards 4-WLT de 
fined for first winglet section only , VIN 
values expected) . 

1-70 

YINUV 

Winglet section inboard surface ordinates 
(VIN values) . 

1-70 

YINLV 

Winglet section outboard surface ordinates 
(VIN values) . 

Cbiit Card set B for 

CASE = 2 (Card 2-A). 

1-10 

BKOD = 1 

. Infinite cylinder (only RADIUS need be 
input) . 


BKOD = - 

1 . Same as BKOD = 1. No embedded body 
grid. Crude grid body representation only. 


BKOD = 2 

!. Simple axisymmetric body definition re- 


quested (input XINB, RIN on card(s) 2-B 
and 3-B) . 



CARD 

CARD 

VARIABLE 


NUMBER 

COLUMN 

NAME 

DESCRIPTION 

Card 1-B 


BKOD = -2 

.Same as BKOD = 2. No embedded body 

(contd) 

11-20 

BKOD = 3. 

BKOD = -3 
BNOSE 

grid. Crude grid body representation only. 
Complex body definition requested (input 
Quick- Geometry model on card(s) 4-B 
through 13- B). 

.Same as BKOD = 3. No embedded body 
grid. Crude grid body representation only. 
Body nose (X- value jl por BKQD = i2 or ±3 


21-30 

BTAIL 

Body tail (X- value) J 


31-40 

BNIN 

Number of axisymmetric body coordinates to 


41-50 

RADIUS 

be input for BKOD = ±2 only (BNIN 1 60). 
Cylinder radius for BKOD = ±1 only. 


51-60 

ANOSB 

ANOSB =0. Sharp nose body^ BKOD = ±2 only 


61-70 

BS 

ANOSB = 1. Blunt nose body J 

Body Cp plot scaling per inch (typically 

NOTE: 

Omit card sets 2-B and 

0.08). 

3-B for BKOD = ±1 or BKOD = ±3. 

Card(s) 2-B 

1-70 

XINB 

Axisymmetric body X coordinates (BNIN 

Card(s) 3-B 

1-70 

RIN 

values) . 

Axisymmetric body radii (BNIN values). 

NOTE: 

Omit card 

sets 4 B through 13 -B for BKOD = ±1 or BKOD = ±2. 

Card 4-B 

1-70 

VTITLE 

Quick- Geometry model title. 

Card 5-B 

1-10 

ACSM 

Number of distinct cross-section models 


(1 ^ ACSM £ 10). (ACSM card sets 6-B 
and 7-B will follow). 
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CARD VARIABLE 


CARD 

NUMBER COLU MN NA ME DESCRI PTION 

Card 6-B 1-10 ADUM Running count of current cross-section 

model (1-ACSM). 

11-20 AARC Number of arcs in current cross-section 

model (1 < AARC < 10). (AARC Card(s) 7-B 
will follow). 

21-60 CTITLE Title or descriptor of current cross-section 

model. 

Card 7-B 1-8 ARCNAM Arc or component name. 

11-14 ASH APE Arc or component shape. 

21-28 PNTNAM(l) Control point name for beginning of this arc. 

31-38 PNTNAM(2) Control point name for termination of this 

arc. 

41-48 PNTMAN(3) Slope control point name for this arc, if 

required. 

Card 8-B 1-10 ANTCSM Number of cross-section models to define 

entire body (1 £ ANTCSM £ 10). (ANTCSM 
card(s) 9-B will follow). 

Card 9-B 1-10 ADUM Running count of current cross-section model 

(1- ANTCSM) 

11-20 AMODEL Index corresponding to already defined 

cross-section models (between 1 and ACSM) . 

21-30 XCSMS1 Starting X-station for current cross-section 

model. 

31-40 XCSMS2 Ending X-station for current cross-section 

model. 

Card 10-B 1-10 BLINE Number of body line models to be defined by 

segments (1 < BLINE £ 25). (BLINE card 
sets 11-B and 12-B follow). 

11-20 ALIAS Number of body line models to be aliased 

(1 1 ALIAS £ 25). (ALIAS card(s) 13-B 
follow) . 
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CARD 

CARD 

VARIABLE 

NUMBER 

COLUMN 

NAME 

description 

Note: 

Card set 11-B and 12 

-B is repeated BLINE times. 

Card 11- B 

1-10 

BLSEG 

Number of segment(s) defining- body line 




model (1 < BLSEG < 10). 


11 

BYORZ 

The letter Y or Z indicates which data 




definition is to follow. 


12-19 

BNAME 

Body line name to be defined. 

Card 12 -B 

1-4 

S SHAPE 

Segment shape. 


11-20 

D(l) 

X-station for beginning of segment. 


21-30 

D(2) 

Y or Z value corresponding to D(l). 


31-40 

D(3) 

X-station for termination of segment. 


41-50 

D(4) 

Y or Z value corresponding to D(3). 


51-60 

D(5) 

X-station for segment slope control point. 


61-70 

D(6) 

Y or Z value corresponding to D(5) 

Note: 

Card set 13 

-B is repeated ALIAS times. 

Card 13 -B 

11 

BYORZ 

The letter Y or Z indicates which data 




definition is to follow. 


12-19 

BNAME 

Body line name to be defined. 


21 

AYORZ 

The letter Y or Z indicates which definition 




is to be used for aliasing. 


22-29 

ANAME 

Body line name to which BNAME is aliased. 

NOTE: 

Read Card i 

set POD, PCD times (Card 3-A). 

Card 1-POD 

1-70 

TITLEP 

Pod/nacelle title to identify graphic and 


printed output. 
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CARD CARD VARIABLE 

NUMBER COL UMN NAME DE S CRIPTIO N 

Card 2-POD 1-10 PTYPE PTYPE = 0. Closed body (pod, tank, etc.)- 

PTYPE = 1. Engine nacelle - cold jet. ) 

j IMFR 

PTYPE =-l. Engine nacelle - hot jet. ) 

and FNPR expected on Card 3-PCD. 

11-20 BNOSE1 Pod nose (X-position) . 

21-30 BTAIL1 Pod tail (X-position) . 

31-40 YBODY1 Pod span position (Y-position) . 

41-50 ZBODY1 Pod height (Z-position). 

51-60 BNIN1 Number of coordinates defining pod 

geometry (BNIN1 ^ 30). 

61-70 ANOSB1 ANOSB1 = 0. Sharp nose pod. 

ANOSB1 = 1. Blunt nose pod. 

Card 3-POD 1-10 BODALF Pod angle -of -attack. 

11-20 BODBET Pod yaw angle (positive-nose of pod away 

from centerline) . 

21-30 FIFP Pod embedded grid code FIFP = 0. Pod 

FIFP = 0. crude grid only. 

FIFP = 1 . Fine pod grid used . 

31-40 FMFR Engine inlet mass flow ratio, MFR. 

41-50 FNPR Engine nozzle pressure ratio, NPR. 

Card 4-POD 1-70 XINBP Axisymmetric pod non-dimensional 

X- ordinates (BNIN1 values). 

Card 5-POD 1-70 RINP Axisymmetric pod non-dimensional radii 

(BNIN1 values). 
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SAMPLE INPUT DATA SETS 


NACA 0012 AIRFOIL 


2, 

0 . 

804 1. 

.94 3, 

.00 

1 . 

56 

1 . 

. 0 0 , 

.25 

0 . 

0 . 

100.0 0 . 



0.0 

0.00200 

0.00500 

0 


0 .05000 

0 . 06000 

0 . 080 00 

0 


0 . 18000 

0.20000 

0.22000 

0 


0 . 32000 

0.34000 

0 . 36000 

0 


0.46000 

0.48000 

0.50000 

0 


0.60000 

0.62000 

0 .64000 

0 


0 .740 00 

0.76000 

0 . 78000 

0 


0.88000 

0.90000 

0 . 92000 

0 


0 . 0 

0.00780 

0.01210 

0 


0 .03540 

0 . 03830 

0 .04300 

0 


0.05610 

0 . 05740 

0 .05840 

0 


0 . 06000 

0 . 05980 

0 , 05940 

0 


0 .05540 

0 . 05420 

0 . 05300 

0 


0 .0457 0 

0 . 04410 

0 . 04230 

0 


0 .03270 

0 .03070 

0 . 02850 

0 


0 .01700 

0 . 01450 

0 .01200 

0 


0 . 0 

-0 .00780 

-0.01210 

-0 


-0.03540 

-0 . 03830 

-0 . 04300 

-0 


-0 .05600 

-0.05730 

-0 .05830 

-0 


-0 . 05990 

-0.05970 

-0 . 05920 

-0 


-0 . 05530 

-0 . 05420 

-0 . 05280 

-0 


-0 .04560 

-0 . 04400 

-0 . 04220 

-0 


-0.03260 

-0 . 03060 

-0 . 02840 

-o 


-0 .01680 

-0 .01440 

-0 .01190 

-o 


150 . 150 . 3 .0 
0 . 1. C .4 
0 . 

01000 0.02000 0.03000 0.04000 
10000 0.12000 0.14000 0.16000 
24000 0,26000 0.2SGC0 0.30000 
38000 0.400CC 0.42000 C. 44 000 
52000 0.54000 0.56000 0.5SCO0 
66000 0.68000 0.7000C 0.729C0 
80000 0.82000 0.84000 0.86000 
94000 0.96000 0.98000 1.00040 
01690 0.02350 0.02830 0.03220 
04670 0.04980 0.05230 0.05440 
05910 0.05960 0.06003 0.06010 
05880 0.05810 0.05730 0.05640 
05170 0.05030 0.04C90 0.04720 
04050 0.03S7C 0.03670 0.03480 
02630 0.02410 0.C2180 0.01940 
00950 0.00680 0.00420 0.00150 
01690 -0.02350 -0.G2830 -0.03220 
04670 -0.04980 -0.05230 -0.05430 
05900 -0.05950 -0.05990 -0.06000 
05860 -0.05S00 -3.05720 -0.35630 
05160 -0.05020 -0.04870 -0.04720 
04040 -0.03860 -0.03660 -0.03470 
02620 -3.02400 -0.02160 -0.01920 
00930 -0.00670 -3.00400 -0.00130 
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NASA AIRFOIL LSI 


2 . 0 . 

722 

- 1.50 4 

. 96 

1 . 38 . 

1.0 0 

. 25 

0 . 0 . 


100.0 0 , 


0 . 0 

0.20000 

0.50000 

1 

7 . 50000 

1 0 . 0000 D 

12.50000 

15 

30.00000 

35.00000 

39.99998 

44 

60.00000 

62.50000 

64 .99998 

67 

77.49998 

80 . 00000 

82.49998 

85 

94.99998 

97.50000 

100.00000 


0.0 

1 .03000 

1.63000 

2 

5.28000 

5.88000 

6 . 37000 

6 

8.20000 

8.40000 

8.49000 

8 

7.67000 

7.39000 

7 . 08000 

6 

4.97000 

4.47000 

3.95000 

3 

1.10000 

0.49000 

-0 . 150 C 0 


0 . 0 

- 0.66000 

- 0 . 97000 

-1 

- 2.94000 

- 3.28000 

- 3.57000 

-3 

- 4.49000 

- 4.52000 

- 4.49000 

-4 

- 3.37000 

- 3. 07000 

- 2 . 760 C 0 

-2 

- I . 10000 

- 0.78000 

- 0 .51000 

-0 

- 0.07000 

- 0.28000 

- 0.71000 



150 . 150 , 3.0 

0 . 1 . 0.4 

0 . 


250 00 
0 0 00 0 
99998 
500 00 
0 0 0 0 0 

2 . 5 0 3 0 C 

17 .50000 
50 .00003 
69. 99998 

87.50000 

3.75000 

19.99998 
55 . 00000 
72.50000 

89.99998 

5 . OCOOO 

25 .00000 
57 .49998 

75.00000 
92.50000 

46000 

3.36000 

4 . 0 C 000 

4.51000 

77000 

46000 

72000 

41000 

7 .12000 
8.33000 
6 . 3 3 0 0 0 
2.85000 

7 .42000 
8 . 07000 
5.91000 
2.28000 

7 .88000 
7.89000 
5.45000 
1 . 7 0 3 0 0 

4 4 0 0 0 
80000 
37000 
430 00 
280 0 C 

- 1.88000 
- 3 .98000 
- 4.17000 
- 2.10000 
- 0.12000 

- 2.23000 
- 4. 150 JO 
- 3.86000 
- 1 . 75000 
0 . 0 

- 2.50000 
- 4 .38000 
- 5.62000 
- 1 . 43 C 3 C 
0.01000 
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KC- 135 

TRANSPORT 

WITH PODS, 

PYLONS, 

3 . 

0.78 

2.0 

5.25 

2. 

1. 

2. 

6 . 

26 . 

0 . 

52 .90 

28 .66 

0 . 0 

70.27 

2.22 

0 . 

0 . 5 

0 . 75 

1 . 25 

10 . 

15 . 

20 . 

25 . 

45. 

50 . 

55. 

60 . 

80 . 

85. 

90 . 

95 . 

0.740 

1.750 

2 . 055 

2.539 

6.572 

7.661 

8.330 

8.860 

8 .766 

8 .258 

7 .620 

6.870 

3.445 

2 . 584 

1 .723 

0.861 

0.000 

-0 .872 

-1 .076 

-1.390 

-3 .665 

-4.510 

-5 . 228 

-5.757 

-6 .044 

-5.678 

-5.185 

-4.618 

-2.309 

-1.732 

-1.155 

-0.577 

32 . 02 

4 .37 

62.89 

2 . 02 

36 . 38 

10.36 

57.86 

1 .75 

0 .536 

1 .399 

1 .654 

2.060 

5.511 

6 .436 

7.070 

7.500 

7.725 

7.382 

6 .895 

6.284 

3 . 188 

2 .391 

1 .594 

0.797 

0.000 

-0.733 

-0 .895 

-1.123 

-2.646 

-3.318 

-3.925 

-4.401 

-5.007 

-4 .793 

-4 .448 

-3 . 995 

-1. 998 

-1.498 

-0.999 

-0.499 

41.41 

16.73 

59.20 

1.46 

0.374 

1 .123 

1 . 340 

1 . 687 

4.699 

5.530 

6 . 047 

6 .422 

6 .736 

6.509 

6.136 

5.643 

2 . 929 

2.197 

1.465 

0 .732 

0.000 

-0.570 

-0 . 684 

-0.837 

-1.787 

-2.265 

-2 .731 

-3.118 

-3 . 782 

-3.696 

-3.487 

-3.174 

-1.593 

-1.194 

-0.796 

-0.398 

50.81 

29.30 

64 . 90 

0.70 

0 .306 

0 . 997 

1.198 

1.519 

4 . 423 

5.210 

5 .704 

6 .050 

6.370 

6 .178 

5.850 

5.412 

2.904 

2. 178 

1 .452 

0 . 726 

0.000 

-0 . 465 

-0.558 

-0.657 

-1.187 

-1.479 

-1 .770 

-2 . 059 

-2.568 

-2 . 501 

-2.330 

-2.078 

-1.039 

-0.775 

-0.520 

-0.260 

69.27 

53.96 

77 , 27 

-1.0 

INBOARD 

2. 

PYLON 
26 . 

1 . 0 

22. 12 

0 . 

-19.0 

1 . 0 

0 . 

0 . 

0 . 5 

0 .75 

1 . 25 

10.0 

15. 0 

20.0 

25 . 0 

45.0 

50 . 0 

55.0 

60.0 

80 . 0 

85 . 0 

90 . 0 

95 . 0 

0 . 

0 . 50 

0 . 90 

1 . 15 

3 .24 

3.86 

4 .34 

4.66 

4.90 

4 .64 

4.30 

3.88 

1 .76 

1 . 22 

0 . 7 

0.26 

0 . 

-0.50 

-0 . 90 

-1.15 

-3.24 

-3.86 

-4.34 

-4 . 66 

-4 . 90 

-4.64 

-4.30 

-3.88 

-1 .76 

-1.22 

-0 . 7 

-0.26 

0 . 

0 . 

1 . 0 

0 . 

OUTBOARD 

2. 

PYLON 
26 . 

1 . 0 

38. 00 

0 . 

-25.0 

1 . 0 

0 . 

0 . 

0 . 5 

0 .75 

1 . 25 

10.0 

15.0 

20 . 0 

25. 0 

45. 0 

50 . 0 

55.0 

60.0 

80 . 0 

85.0 

90.0 

95.0 

0 . 

0 .50 

0 . 90 

1 . 15 

3.24 

3.86 

4 . 34 

4 . 66 

4 .90 

4.64 

4.30 

3 . 88 

1 . 76 

1.22 

0.7 

0 . 26 

0 . 

-0.50 

-0 . 90 

-1.15 

-3.24 

-3.86 

-4.34 

-4 .66 

-4.90 

-4 .64 

-4 .30 

-3.88 

-1.76 

-1 . 22 

-0 . 7 

-0.26 

0 . 

0 . 

1 . 0 

0 . 


AND WINGLETS 


100 . 

80 . 

-4.25 
1 . 0 

1440 . 0 

2.50 

5.0 

30 . 

35 . 

65 . 
100 . 

70 . 

3.470 

4.830 

9.137 

9.218 

6 . 029 
0 . 

5.168 

-1 . 915 

-2.613 

-6.104 

-6.275 

-4.041 
0 . 

0. 

1 . 

-3.464 

2.840 

3.998 

7.760 

7 . 920 

5.570 
0 . 

4.782 

-1.475 

-1 . 915 

-4 .746 

-4 . 980 

-3.496 
0 . 

1 . 0 

-2. 996 

2.350 

3.360 

6.666 

6.805 

5.049 
0 . 

4.380 

-1 .053 

-1 . 307 

-3 . 425 

-3.642 

-2.787 
0 . 

1 . 0 

-2 . 389 

2.160 

3. 130 

6.280 

6.430 

4 .880 
0. 

4.267 

-0.750 

-0.895 

-2.275 

-2.445 

-1.818 

0. 

0. 

-1.558 


1 . 0 


2 .50 

5.0 

30 . O 

35 . 0 

65.0 
100 . 0 

70 . 0 

1.70 

2.36 

4. 90 

5.02 

3.40 

0 . 

2.88 

-1.70 

-2.36 

-4.90 

-5.02 

-3.40 

0 . 

0 . 

1 . 0 

-2.88 

2.50 

5.0 

30 . 0 

35 . 0 

65. 0 

100.0 

70.0 

1.70 

2.36 

4.90 

5 .02 

3 .40 
0 . 

2.88 

-1.70 

-2.36 

-4.90 

-5 . 02 

-3.40 

0 . 

0 . 

-2.88 


1 . 

0.8 

7.5 
40 . 

75. 

5.813 

9.105 

4.300 

“3.178 

-6.250 

- 2.886 


4 .850 
7.910 
3 . 985 

-2.290 
-5.07 0 
-2.497 


4.114 

6.829 

3.687 

-1 . 544 
-3.763 
-1.991 


3.857 

6.449 

3.605 

-1.041 
-2.543 
-1 .299 


7 . 5 
40.0 
75 . 0 

2.82 

5.02 

2.34 

-2.82 
-5 . 02 
-2.34 


7.5 
40 . 0 
75 . 0 

2.82 
5. 02 
2.34 

- 2.82 
-5.02 
-2 .34 
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KC- 135 UINGLET 


2. 

38. 

1 . 0 

75 . 0 




72.07 

53 . 96 

77.27 

-9 . 0 

1 . 0 



0 . 

0 . 2 

0.5 

1 . 25 

2 .50 

3.75 

5.00 

7 .50 

10.0 

12.5 

15.0 

17 . 5 

20.0 

25.0 

30 . 0 

35.0 

90.0 

95 . 0 

50 . 0 

55 . 0 

57 . 5 

60.0 

62 . 5 

65 . 0 

67.5 

70.0 

72 . 5 

75.0 

77 .5 

80.0 

82.5 

85.0 

87 . 5 

90 . 0 

92 . 5 

95 . 0 

97 . 5 

100.0 





0 . 

0.77 

1.19 

1.79 

2.99 

2 . 96 

3.33 

3.89 

9 . 33 

9.69 

9 . 99 

5.25 

5.47 

5.81 

6.05 

6 .21 

6 . 28 

6 .27 

6 .18 

3. 99 

5.87 

5 . 72 

5.59 

5.33 

5.08 

9 .81 

4.51 

4.19 

3.89 

3 .99 

3. 11 

2.70 

2.28 

: .84 

1 . 38 

0.89 

0 .38 

-0 .20 





0 . 

-0 . 31 

-0.91 

-0.60 

-0.77 

-0.90 

-1.00 

-1.18 

-1.32 

-1.99 

-1.59 

-1.61 

-1.67 

-1.75 

-1.76 

-1.79 

-1.68 

-1.55 

-1.99 

-1.22 

-1.06 

-0 . 90 

-0.71 

-0.52 

-0.33 

-0 . 15 

0 . 04 

C .20 

0 . 36 

0.99 

0.60 

0 .65 

G . 64 

0 .59 

0 . 45 

0 .21 

-0.13 

-0 .67 





78 . 32 

61.96 

80.00 

-9 . 0 

0 . C 



-2 . 

0 . 00 

108.30 

22. 

5 . 09 

l 4 

. 4 

0 . 

2 . 0 

9 . 0 

6 . 0 

8.0 

10 . 0 

12.0 

19.0 

16 . 0 

20 . 0 

30 . 0 

90 . 0 

50 . 0 

60.0 

70.0 

78 . 0 

80.0 

85 . 0 

90 . 0 

95 . 0 

100.0 

1C8.3 







0 . 

2 .50 

3 . 35 

3.90 

9 .30 

4 .60 

4.80 

9 . 95 

5 . 09 

5 . 09 

5.09 

5.09 

5 .04 

5.04 

5 . 09 

5 . 09 

9 . 92 

9.65 

4 .20 

3 .55 

2 . 65 

2.60 







KC 1 35 A 

INBOARD ENGINE NACELLE 





0 . 

33 . 969 

95.610 

22.05 

-3.816 

17 . 

0 . 

2 . 

0 . 

0 . 

0 . 

0 . 



0 . 

0.58 

1.81 

9.61 

7 . 57 

10.37 

16.14 

21.90 

27 .66 

39.19 

93.31 

50.72 

6 2.24 

79.53 

91.06 

96.82 

100.0 





9.72 

10.21 

10.87 

11.69 

12.35 

12.84 

13.50 

13.83 

19 .16 

19 .00 

19.16 

14.00 

13.67 

11 . 94 

10.09 

8.90 

7 .91 





KC 1 35 A 

OUTBOARD 

ENGINE NACELLE 




0 . 

95.890 

57 . 986 

38 . 15 

-3.816 

17 . 

0 . 

2 . 

0 . 

0 . 

0 . 

0 . 



0 . 

0 . 58 

1.81 

9.61 

7 . 57 

10.37 

16.14 

21.90 

27.66 

39. 19 

93.31 

50.72 

62 . 24 

79.53 

91.06 

96.82 

100.0 





9 . 72 

10.21 

10.87 

11.69 

12.35 

12.84 

13.50 

13.83 

19 . 16 

19.00 

19 . 16 

14.00 

13.67 

11.94 

10.09 

8.90 

7.91 
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GULFSTREAM III CONFIGURATION WITH WINGLET AND NACELLE 


3 . 0 . 

780 4 , 

.00 3 . 

, 10 

100 . 

80 . 

1.0 

U . I . 

14 . 26 

X ■ 
1 . 

!o 451.3 

- 47.0 

134582.4 

0.6 

275.0 0 . 

508.72 0 , 

,0 

1.0 



0 . 0 

0 .00500 

0.00750 

0 .01250 

0.02500 

0.05000 

0 . 07500 

0.10000 

Q '15000 

0.20000 

0.25000 

0.30000 

0.35000 

0 . 4 U 000 

0.45000 

0.50000 

0.55000 

0.60000 

0.65000 

0.70000 

0 .75000 

0.80000 

0.85000 

0.90000 

0 . 95000 

1.00000 



- 0.05467 

- 0 . 04646 

- 0 . 04435 

-0 . 04088 

- 0.03409 

- 0.02419 

- 0 .01713 

- 0.01245 

- 0.00905 

- 0.00845 

- 0.00902 

- 0.01059 

- 0.01369 

- 0 .01760 

- 0.02274 

- 0.02857 

- 0.03483 

-0 . 04147 

- 0.04846 

- 0 .05585 

-0 . 06385 

- 0.07309 

- 0.08260 

- 0 .09212 

- 0.10163 

- 0.11114 



- 0.05467 

- 0.06342 

- 0.06548 

- 0.06876 

- 0.07480 

- 0.08334 

- 0.08974 

- 0.09489 

- 0.10287 

- 0.10912 

- 0.11465 

- 0.11965 

- 0.12419 

- 0.12819 

- 0.13145 - 

- 0.13344 

- 0 .13438 

- 0.13447 

- 0.13398 

- 0.13299 

- 0 .13151 

- 0.12940 

- 0.12634 

- 0.12234 

- 0.11750 

- 0.11150 



302.7 45 

.0 519.7 0 . 

, 0 

1 . 0 



- 0 . 06251 

- 0.05440 

-0 . 05240 

- 0 . 04908 

- 0.04276 

- 0.03373 

-0 . 02726 

- 0 . 02281 

- 0.01835 

- 0 .01672 

- 0.01651 

- 0.01725 

- 0.01905 

-0 . 02224 

- 0 . 02643 

- 0 .03161 

- 0.03744 

- 0 .04390 

-0 . 05090 

- 0.05348 

-0 . 06670 

- 0.07555 

- 0.08430 

- 0 .09277 

- 0 .10112 

- 0.10947 



- 0. 06251 

- 0.07061 

-0 . 07249 

- 0.07551 

- 0 . 08106 

-0 . 08880 

- 0 .09450 

-0 . 09900 

- 0.10588 

- 0 .11137 

- 0.11608 

- 0.12020 

- 0.12382 

- 0.12697 

- 0.12963 

- 0.13143 

- 0.13217 

- 0 .13214 

- 0.13149 

- 0.13026 

- 0 .12848 

- 0 . 12603 

- 0 . 12275 

- 0.11895 

- 0.11473 

- 0.10984 



321.2 75 

.0 526.9 0 . 

. 0 

1.0 



- 0.06846 

- 0.06048 

- 0.05856 

- 0 .05540 

- 0.04943 

- 0.04099 

- 0.03498 

- 0.03067 

- 0.02556 

- 0.02323 

- 0,02245 

- 0.02275 

- 0.02382 

-0 . 026 18 

- 0.02970 

- 0.03434 

-0 .03976 

- 0.04588 

- 0.05279 

- 0.06044 

-0 . 06884 

- 0.07764 

- 0.08605 

- 0 .09375 

- 0.10098 

- 0.10821 



- 0.06846 

- 0.07600 

-0 . 07775 

-0 . 08055 

- 0 . 08568 

- 0.09275 

- 0.09783 

- 0 .10172 

- 0.10771 

- 0 .11242 

- 0 .11640 

- 0.11986 

- 0.12294 

- 0.12573 

- 0.12822 

- 0.12984 

- 0 .13047 

- 0 .13033 

- 0.12952 

- 0.12813 

-0 . 12612 

- 0.12339 

- 0.12015 

- 0 .11659 

- 0.11277 

- 0.10858 



328.0 86 

.0 529.6 0 . 

, 0 

1.0 



- 0.07081 

- 0 .06289 

- 0.06099 

- 0.05790 

- 0.05208 

- 0.04391 

- 0.03804 

- 0.03377 

- 0.02847 

-0 . 02586 

- 0.02487 

- 0 . 02497 

- 0.02590 

-0 . 02787 

- 0.03113 

- 0.03550 

- 0 . 04072 

- 0 ,04669 

- 0.05350 

- 0.06110 

-0 . 06947 

- 0 . 07817 

- 0.08637 

- 0 .09368 

- 0.10070 

- 0.10771 



- 0.07080 

- 0.07812 

- 0.07981 

- 0.08250 

- 0.08746 

- 0.09424 

- 0.09900 

- 0.10266 

- 0.10825 

- 0.11263 

- 0 .11633 

- 0.11957 

- 0.12250 

-0 . 12519 

- 0.12767 

- 0.12922 

- 0.12980 

- 0.12961 

- 0.12874 

- 0.12727 

- 0.12515 

- 0.12235 

- 0.11915 

- 0.11571 

- 0.11204 

- 0.10808 



352.7 126.0 539 . A 0 , 

.0 

1.0 



- 0.08021 

- 0.07258 

-0 . 07083 

- 0.06805 

- 0.06282 

- 0.05561 

- 0.05032 

- 0 . 04618 

- 0.04025 

- 0.03671 

- 0 .03486 

- 0.03425 

- 0.03455 

- 0.03554 

- 0.03751 

- 0.04070 

- 0.04509 

- 0.05026 

- 0.05626 

- 0.06320 

-0 . 07088 

- 0.07824 

- 0 . 08511 

- 0.09198 

- 0.09885 

- 0.10572 



- 0.08021 

-0 . 08648 

- 0.08793 

- 0.09022 

- 0.09433 

- 0.09957 

- 0.10288 

- 0.10532 

- 0 .10913 

- 0.11228 

- 0 .11516 

- 0.11789 

- 0.12051 

- 0 .12308 

- 0 . 12539 

- 0.12671 

- 0 .12711 

- 0.12666 

- 0 , 12551 

- 0.12366 

- 0.12116 

- 0.11827 

- 0 .11531 

- 0 .11229 

- 0 .10922 

- 0 .10608 



364.4 145.0 544.0 0 

.0 

1.0 



. - 0.08522 

- 0.07780 

- 0.07614 

- 0 .07350 

- 0.06866 

- 0.06191 

- 0.05687 

-0 . 05281 

- 0.04666 

- 0.04270 

- 0.04042 

- 0.03942 

- 0.03938 

- 0 .04008 

- 0.04145 

- 0.04397 

- 0.04775 

- 0 .05246 

- 0.05778 

- 0.06389 

- 0 ,07057 

- 0.07739 

- 0.08420 

- 0.09102 

- 0 .09784 

- 0.10465 



-0 . 08522 

- 0 .09085 

- 0.09215 

- 0.09418 

- 0.09766 

- 0.10159 

- 0 .10387 

- 0 .10566 

- 0.10872 

- 0.11153 

- 0.11420 

- 0.11679 

- 0 .11938 

-0 . 12194 

-0 . 12415 

- 0.12533 

-0 . 12563 

- 0.12505 

-0 . 12369 

- 0.12162 

- 0.11899 

- 0.11620 

- 0.11340 

- 0.11061 

- 0.10781 

-0 . 10502 



378.0 167.0 549.3 0 

. 0 

1.0 



- 0.08537 

- 0.07765 

-0 . 07592 

- 0.07324 

- 0.06830 

- 0.06147 

-0 . 05641 

- 0.05232 

- 0.04615 

- 0.04216 

- 0.03986 

- 0.03883 

- 0 .03878 

-0 . 03948 

-0 . 04089 

- 0 . 04336 

-0 . 04701 

- 0.05165 

- 0 . 05690 

-0 . 06290 

- 0.06947 

- 0.07624 

- 0.08301 

-0 . 08978 

- 0.09655 

- 0.10333 



- 0.08537 

- 0.09119 

- 0.09249 

- 0.09450 

- 0.09783 

- 0.10149 

- 0.10387 

- 0 .10575 

- 0.10888 

- 0.11169 

- 0.11433 

- 0.11687 

- 0.11939 

- 0.12186 

-0 . 12382 

- 0.12480 

- 0.12496 

- 0.12426 

- 0.12281 

- 0.12070 

- 0 .11797 

- 0.11511 

- 0.11226 

- 0.10940 

- 0.10654 

- 0.10369 



402.7 207.0 559.0 0 

. 0 

1.0 



- 0.08567 

- 0.07738 

- 0.07555 

- 0.07270 

- 0.06756 

- 0.06057 

- 0.05544 

- 0.05130 

- 0.04508 

- 0.04104 

- 0 .03868 

- 0.03759 

- 0 .03753 

-0 . 03824 

- 0.03970 

- 0.04208 

-0 . 04556 

- 0.05004 

- 0.05508 

- 0 . 06036 

- 0. 06721 

- 0.07387 

- 0.08054 

- 0.08721 

- 0.09388 

- 0 .10055 



- 0.08567 

- 0.09187 

-0 . 0931 9 

- 0 .09518 

- 0.09825 

- 0.10147 

- 0.10391 

- 0 .10592 

- 0.10919 

- 0.11198 

- 0.11453 

- 0.11691 

- 0.11928 

- 0.12150 

-0 . 12294 

- 0 . 12361 

-0 . 12349 

- 0.12263 

-0 . 12099 

- 0.11872 

- 0.11582 

- 0 .11283 

- 0.10985 

- 0 .10687 

- 0 .10389 

- 0.10091 
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452.7 

288.0 

578.7 0 

. 0 

1.0 



-0 . 08651 

- 0 . 07666 

- 0.07457 

- 0.07130 

- 0.06566 

- 0.05815 

- 0 . 05279 

- 0.04849 

- 0 . 04213 

- 0 .03799 

- 0 .03548 

- 0 .03425 

- 0 .03409 

- 0 . 03478 

- 0.03643 

- 0 . 03865 

- 0.04176 

- 0 . 04568 

- 0.05039 

- 0.05555 

- 0 . 06134 

-0 . 06751 

- 0 . 07386 

-0 . 08021 

- 0 . 08656 

- 0.09291 



-0 . 08651 

- 0 .09375 

- 0 .09513 

- 0.09705 

- 0.09957 

- 0.10218 

- 0.10442 

- 0.10635 

- 0.10954 

- 0.11212 

- 0.11434 

- 0.11631 

- 0.11809 

- 0.11931 

- 0.11993 

- 0.11992 

- 0 .11917 

- 0 .11784 

- 0.11582 

- 0.11319 

- 0.10984 

-0 . 10653 

- 0.10321 

- 0 .09990 

- 0.09659 

- 0 .09327 



479.9 

332.0 

589.4 0 . 

. 0 

1.0 



-0 . 08717 

- 0 . 07612 

-0 . 07384 

- 0 . 07037 

- 0.06425 

- 0.05630 

- 0.05073 

- 0.04631 

-0 . 03987 

- 0.03561 

- 0.03299 

- 0.03162 

- 0.03139 

- 0 . 03213 

- 0.03389 

- 0 .03608 

- 0.03890 

- 0.04252 

- 0.04690 

- 0.05181 

- 0.05713 

-0 . 06280 

- 0.06885 

-0 . 07489 

-0 . 08094 

- 0.08698 



-0 . 08717 

- 0 .09518 

- 0.09664 

- 0.09860 

- 0 .10081 

- 0.10313 

- 0.10510 

- 0 .10676 

- 0 .10946 

- 0.11164 

- 0.11345 

- 0 .11508 

- 0.11625 

- 0. 11699 

- 0.11714 

- 0.11670 

- 0 .11564 

- 0.11393 

- 0.11167 

- 0.10869 

- 0 .10516 

- 0.10160 

- 0 .09803 

- 0.09447 

- 0.09091 

- 0.08734 



502.8 

369.0 ! 

598.4 0 . 

0 

1.0 



- 0.08789 

- 0 .07569 

- 0.07323 

- 0.06933 

-0 . 06280 

- 0.05435 

- 0.04854 

-0 . 04396 

- 0 .03737 

- 0.03300 

- 0.03022 

- 0.02875 

- 0.02840 

- 0 . 02925 

- 0 . 03101 

-0 . 03316 

- 0.03587 

- 0 .03923 

- 0.04321 

- 0.04787 

-0 . 05279 

-0 . 05799 

- 0.06352 

-0 . 06915 

-0 . 07479 

- 0.08042 



- 0.08789 

- 0.09680 

- 0.09832 

- 0.10030 

- 0.10231 

- 0.10410 

- 0 .10551 

- 0.10674 

- 0.10877 

- 0.11040 

- 0.11186 

- 0.11299 

- 0.11374 

- 0.11403 

- 0.11375 

- 0.11291 

- 0.11147 

- 0.10950 

- 0.10693 

- 0.10364 

- 0.09992 

- 0.09609 

- 0.09227 

- 0.08844 

- 0.08461 

- 0.08078 



513.9 

387.0 602.8 0 . 

0 

1.0 



- 0.08833 

- 0.07543 

- 0.07286 

-0 . 06881 

- 0.06191 

- 0.05316 

- 0 . 04718 

- 0 . 04257 

- 0.03587 

- 0.03146 

- 0.02861 

- 0.02701 

- 0.02658 

- 0.02753 

- 0 . 02929 

-0 . 03145 

- 0 . 03408 

- 0.03728 

- 0.04109 

- 0.04556 

- 0.05032 

- 0.05533 

- 0.06049 

- 0.06582 

- 0.07115 

- 0 .07648 



-0 . 08833 

- 0.09769 

- 0.09931 

- 0.10130 

- 0.10324 

- 0.10445 

- 0.10545 

-0 . 10638 

- 0 .10792 

-0 . 10937 

- 0.11060 

- 0.11153 

- 0.11204 

- 0 .11203 

- 0.11154 

- 0.11053 

- 0.10894 

- 0.10678 

- 0.10396 

- 0.10061 

- 0 .09676 

- 0.09278 

-0 . 08879 

- 0.08481 

- 0.08083 

- 0.07684 



540.5 430.0 613.2 0 . 

0 

1,0 



-0 . 08969 

-0 . 07483 

- 0 .07178 

- 0.06722 

- 0.05948 

- 0.04968 

- 0.04313 

- 0.03819 

- 0.03115 

- 0 .02667 

- 0.02347 

- 0.02156 

- 0.02091 

- 0 . 02191 

-0 . 02369 

- 0.02596 

-0 . 02852 

- 0.03143 

- 0,03481 

- 0.03863 

-0 . 04305 

- 0 . 04757 

- 0.05200 

-0 . 05614 

- 0.06014 

- 0.06414 



- 0.08969 

- 0.10053 

- 0.10233 

- 0 .10458 

- 0 . 10632 

- 0.10538 

- 0 .10458 

- 0.10425 

- 0 .10443 

-0 . 10528 

- 0.10594 

- 0.10620 

- 0.10601 

- 0 .10540 

- 0 .10432 

- 0.10266 

- 0.10052 

- 0.09781 

- 0.09462 

- 0.09110 

-0 . 08673 

- 0 . 08229 

- 0.07784 

- 0.07339 

- 0.06894 

- 0.06450 



552.8 450.0 

618.1 0 

.0 

1.0 



- 0.09055 

- 0.07451 

- 0.07127 

- 0 . 06635 

- 0.05796 

- 0.04758 

- 0.04061 

- 0.03544 

- 0.02824 

-0 . 02363 

- 0.02027 

- 0.01803 

- 0.01754 

-0 . 01840 

- 0.02010 

-0 . 02231 

- 0.02490 

- 0 .02782 

- 0.03104 

- 0.03459 

- 0.03870 

- 0.04311 

- 0.04732 

- 0.05111 

- 0.05371 

- 0.05632 



- 0 .09055 

- 0.10234 

- 0.10428 

- 0.10664 

- 0.10846 

- 0.10626 

- 0 .10401 

- 0. 10236 

-0 .10161 

- 0 . 10218 

- 0 .10245 

- 0.10231 

-0 . 10179 

- 0 .10087 

- 0.09939 

- 0.09748 

- 0.09503 

- 0.09190 

- 0.08859 

- 0.08503 

- 0 . 08030 

- 0 .07558 

-0 . 07085 

-0 . 06613 

- 0 . 06140 

- 0.05668 




G- III WINGLET 


3 . 

24 . 

1 . 0 

75.0 




576.0 

450.0 

618.1 

0 . 0 

1.0 



0 . 

0 .005 

0.01 

0 . 02 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 

0 . 95 

1 . 0 





0.0 

0.0088 

0.01198 

0 .01650 

0.02580 

0.03615 

0.04311 

0.04762 

0 . 05049 

0 .05175 

0 . 05150 

0.05020 

0 . 04805 

0 ,04534 

0.04212 

0.03850 

0 .03425 

0.02990 

0.02515 

0.02034 

0 .01535 

0.01060 

0 .00605 

0.00170 





0.0 

- 0. 00780 

- 0.01007 

- 0.01290 

- 0.01740 

- 0.02163 

- 0 . 02424 

- 0.02577 

- 0 . 02650 

- 0.02670 

- 0 . 02640 

- 0.02555 

- C . 02440 

- 0 . 02295 

- 0.02100 

- 0.01845 

- 0.01500 

- 0.01070 

- 0.00620 

- 0 .00250 

0 . 00010 

0 .00155 

0.00100 

- 0.00220 





602.5 

483.5 

631 . 0 

0 . 

1.0 



0 . 0 

0 .01166 

0 .01608 

0.02183 

0.03160 

0 . 04160 

0 .04790 

0.05270 

0 .05625 

0 . 05855 

0.05995 

0.06060 

0.06060 

0 . 05967 

0 . 05790 

0 ,05540 

0 .65200 

0 . 04755 

0.04230 

0.03665 

0.03030 

0.02300 

0.01385 

0 .00300 





0.0 

- 0.00757 

- 0.01012 

- 0.01300 

- 0.01709 

- 0 . 02045 

- 0 .02202 

- 0.02265 

- 0.02265 

- 0.02220 

- 0.02115 

- 0.01975 

- 0.01793 

- 0.01555 

- 0.01280 

- 0 .00973 

- 0.00630 

- 0.00270 

0.00070 

0.00360 

0 .00480 

0.00430 

0.00200 

- 0.00220 





629.0 

517 . 0 

644.0 

0.0 

0.0 



- 3 . 

8.0 

864.5 

0 . 

0 . 

0 . 

0.4 
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GULFSTREAM FUSELAGE QUICK-GEOMETRY MODEL 
3 . 


1 . 

2. 

NOSE TO 

WINDSHIELD 

BASE 



BODYLO 

ELLI 

BDYBCL 

BDYMHB 

BDYLSCP 



BODYHI 

ELLI 

BDYMHB 

BDYTCL 

BDYUSCP 



2. 

3. 

WINDSHIELD 




CANLO 

ELLI 

BDYBCL 

BDYMHB 

BDYLSCP 



CANSI 

ELLI 

BDYMHB 

CANLOW 

CANLSCP 



WINDF 

LINE 

CANLOW 

BDYTCL 




3. 

4. 

CANOPY 





CANOPLO 

ELLI 

BDYBCL 

BDYMHB 

BDYLSCP 



CANOPSI 

ELLI 

BDYMHB 

CANLOW 

CANLSCP 



WINDSI 

LINE 

CANLOW 

CANHIE 




CANOrUP 

ELLI 

CANHIE 

BDYTCL 

CANTSCP 



5 . 
1 . 

1. 

8.0 

64 . 0 




2. 

2. 

64.0 

83.5 




3. 

3. 

83.5 

133.0 




A. 

1 . 

133.0 

575.5 




5. 

1 . 

575.5 

864.5 




10 . 

11 . 






3. 

2BDYBCL 




8.0 

53. 0 

El LX 

8.0 

78.5 

133.0 

53.0 

LINE 

133 . 0 

53 . 0 

575.5 

53.0 


53.0 

CUBI 

575.5 

53.0 

864.5 

114.5 

625.0 

5 . 

ZBDYTCL 





94.5 

EL LX 

8.0 

78.5 

64.0 

113.5 

8.0 

LINE 

64.0 

113.5 

83.5 

132.5 


147 . 0 

CUBI 

83.5 

132.5 

133.0 

147.0 

98.0 

LINE 

133 . 0 

147 . 0 

575 . 5 

147 . 0 


147 . 0 

CUBI 

575.5 

147.0 

864.5 

114.5 

658.0 

I . 

YCENTER 






LINE 

8.0 

0.0 

864.5 

0 . 0 



3 . 

ZBDYMHB 




135.5 

100.0 

CUBI 

8 . 0 

78.5 

181 . 0 

100.0 

LINE 

181 . 0 

100.0 

502.0 

100.0 


100 . 0 

CU3 1 

502.0 

100.0 

864.5 

114.5 

616.5 

3 . 

YBDYMHB 





47 . 0 

EL LX 

8.0 

0 . 

206 . 0 

47 . 0 

8.0 

LINE 

206 . 0 

47 . 0 

502.0 

47 . 0 


47 . 0 

CUBI 

502.0 

47 . 0 

864.5 

0. 

636.0 

1 . 

ZCANLOW 






LINE 

64 . 0 

113.5 

133.0 

113.5 



I . 

ZCANHIE 






LINE 

83.5 

132.5 

133.0 

132.5 



2. 

YCANLOW 






LINE 

64 . 0 

0 . 

79.5 

28 . 0 



LINE 

79.5 

28.0 

133.0 

42.0 



2. 

YCANHIE 






LINE 

83.5 

0 . 

96 . 0 

24.5 



LINE 

96 . 0 

24.5 

133 . 0 

34 . 0 



1 . 

YCANTSCP 






LINE 

83.5 

0 . 

133.0 

22.0 




ZMAPAXIS 

ZBDYMHB 






YMAP AXIS 

YCENTER 






ZCANTSCP 

ZBDYTCL 






ZCANLSCP 

ZCANLOW 






YCANL5CP 

YBDYMHB 






YBDYTCL 

YCENTER 






YBDYBCL 

YCENTER 






YBDYLSCP 

YBDYMHB 






ZBDYLSCP 

ZBDYBCL 






YBDYUSCP 

YBDYMHB 






ZBDYUSCP 

ZBDYTCL 





G-III NACELLE 




16.0 


1 . 

481.25 

668.75 

78.95 

73.3 

0 . 0 

2.75 

-0.5 

0 . 

0.66 

2.0 



0 . 

6.67 

13.33 

20 .00 

26.67 

33.33 

40.00 

46.67 

53.33 

60.00 

66.67 

73.33 

80.00 

86.67 

93.33 

100 . 0 






9.82 

12.0 

13.05 

13.69 

13.89 

14.04 

14 . 04 

14.04 

13.96 

13.62 

13.13 

12.21 

11.16 

9.89 

8.42 

6.81 
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C-5A TRANSPORT WITH 

ENGINE 

PODS AND 

PYLONS 



3 . 

0.775 

2.50 

2.5 

100 . 

80 . 

3.0 

2 . 

0 . 

2. 





5. 

31 . 

1.0 

165.3 

20 . 0 

18720 . 0 

1 . 0 

109.96 

0 . 

189.17 

4.18 

0 . 




0 . 0 

0 . 25000 

0 .5000 0 

0.75000 

1.00000 

2.00000 

4.00 00 0 

6.00000 

8.00000 

10.00000 

12.00000 

14. 0 0 0 0 0 

16.00000 

18.00000 

20.00000 

25.00000 

30.00000 

35 .00000 

40.00000 

45.00000 

50.00000 

55.00000 

60.00000 

65.00000 

70 . 00000 

75.00000 

80 .00000 

85.00000 

90.00000 

95.00000 

100 . 0000Q 





0 . 56298 

1 . 16181 

1.49063 

1.74381 

1 . 95621 

2.59646 

3.354 04 

3 . 94282 

4.39117 

4.76414 

5 .05056 

5.33434 

5.58457 

5.79552 

6 .00648 

6 .39246 

6.67431 

6.85000 

6 . 91854 

6 .91625 

6.76694 

6 . 55966 

6 .26794 

5 . 90356 

5.45495 

4. 90702 

4.25141 

3.48790 

2.56929 

1 .46900 

0 . 12887 





0 . 56336 

-0.41909 

-0.70382 

-0 .89820 

-1.06155 

-1 . 56102 

-2.31165 

-2 . 92237 

-3.41914 

-3 .83526 

-4 .15688 

-4.47850 

-4.76203 

-5.00587 

-5.24884 

-5.68123 

-5 . 98528 

-6 .16150 

-6.20322 

-6 .11709 

-5.90316 

-5.57466 

-5.13838 

-4.58723 

-3 . 96548 

-3.29717 

-2.57890 

-1.84615 

-1 . 16210 

-0 . 57491 

-0 .12896 





147.37 69 

.5 201.04 2, 

o 

o 

. 0 



0 . 58099 

1 . 16726 

1 .49248 

1.75712 

1 . 96323 

2.60530 

3 . 38237 

3.98024 

4.44052 

4.80000 

5.11966 

5 .41050 

5.66225 

5.87825 

6 .09324 

6.49263 

6.78225 

6 . 97204 

7 .06200 

7 .06200 

6 . 97212 

6 . 78242 

6 .50284 

6 . 14336 

5 . 68417 

5. 11525 

4.41675 

3.53518 

2.53115 

1 . 44331 

0.11019 





0 . 58099 

-0 .39277 

-0.70676 

-0 . 90126 

-1 .06083 

-1 .51522 

-2.11367 

-2.57321 

-2.92607 

-3.20495 

-3.40705 

-3.60305 

-3.77372 

-3.90973 

-4 . 04572 

-4.2S3S8 

-4 . 42525 

-4 .49639 

-4 . 48732 

-4 .40760 

-4.24770 

-4.01791 

-3.71824 

-3 . 34869 

-2.90013 

-2.42997 

-1 . 94043 

-1.46043 

-0.98114 

-0 . 53175 

-0.11019 





154.10 82 

. 0 2 

03.17 2. 

.68 1 

. 0 



0 . 53031 

1 . 15696 

1 .49673 

1.75361 

1 . 96380 

2.60772 

3.38790 

3.98975 

4.43000 

4.79000 

5. 14 0 97 

5 .43662 

5.69037 

5.90815 

6 . 12201 

6 .53237 

6.82420 

7.00831 

7.10630 

7 .11697 

7 . 028 03 

6 .34322 

6.57149 

6.20261 

5 .75804 

5 . 13479 

4 . 47149 

3.59955 

2.59158 

1 . 37545 

0.10949 





0 . 58031 

-0.39956 

-0.70439 

-0.89548 

-1.05719 

-1 .49729 

-2.06062 

-2 . 47998 

-2 .79766 

-3 . 03450 

-3.20384 

-3 . 36684 

-3 . 50283 

-3.60234 

-3.70284 

-3.88321 

-3 . 98274 

-4 .01792 

-3 . 996 96 

-3 . 92039 

-3.77952 

-3.56939 

-3.31939 

-2 . 98628 

-2 .60584 

-2 . 19457 

-1.77261 

-1 . 35957 

-0 . 94637 

-0 . 52253 

-0 . 10949 





167.59 108.0 211.27 2, 

.18 1 

. 0 



0.60323 

1 .24885 

1 . 57273 

1.82851 

2.04443 

2.69783 

3.47276 

4 . 06878 

4.51000 

4.89000 

5.21000 

5.51035 

5.74338 

5.96338 

6.17917 

6 . 57927 

6.85882 

7 . 04174 

7.11491 

7 .10258 

6 . 98062 

6.76147 

6 .45772 

6 .05988 

5.58484 

5.00134 

4.28537 

3.44006 

2.47901 

1 . 35890 

0.11050 





0 .60323 

-0 . 38175 

-0 . 72092 

-0.91419 

-1 .07370 

-1 .52129 

-2.06233 

-2.45674 

-2.74150 

-2.95825 

-3. 10856 

-3.25177 
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C-141 TRANSPORT WITH PODS AND 
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LOCKHEED 

C-Hl QUICK-GEOMETRY FUSELAGE 

MODEL 



4 . 
1 . 

2. 

NOSE TO TAIL 




BODYLO 

ELLI 

BDYBCL 

BDYNHB 

BDYLSCP 



BODYHI 

ELLI 

BDYNHB 

BDYTCL 

BDYUSCP 



2 . 

3 . 

FUSELAGE 

WITH WING 

MOUNT FAIRING 


FAIRLO 

ELLI 

BDYBCL 

BDYNHB 

BDYLSCP 



FAIR5I 

LINE 

BDYNHB 

BS 




FAIRTP 

LINE 

BS 

BDYTCL 




3 . 

5. 

FUSELAGE 

WITH WHEEL AND WING 

FAIRING 


WHEELO 

LINE 

BDYBCL 

WL 




WHEEL 51 

LIME 

WL 

WT 




WHEELUP 

LIME 

WT 

BDYNHB 




FAIRS 

LINE 

BDYNHB 

BS 




FAIRT 

LINE 

BS 

BDYTCL 




4 . 

4 . 

FUSELAGE 

WITH WHEEL FAIRING 



WHELO 

LINE 

BDYBCL 

WL 




WHELS 

LINE 

WL 

WT 




WHELU 

LINE 

WT 

BDYNHB 




UPPER 

ELLI 

BDYNHB 

BDYTCL 

BDYUSCP 



5 . 
1 . 

1 . 

230.4 

650 . 0 




2 . 

2. 

650 . 0 

794.0 




3 . 

3 . 

794.0 

1 050.0 




4 . 

4 . 

1050.0 

1198.0 




5 . 

1 . 

1198.0 

1822 . 0 




11 . 

8 . 






3 . 

ZBDYBCL 






EL LX 

230.4 

-15.0 

445.0 

-85 . 0 

230.4 

-85 . 0 

LINE 

445.0 

-85.0 

1100.0 

-85. 0 



CUBI 

1100 . 0 

-85.0 

1822.0 

70.0 

1485 . 0 

-85 . 0 


I 
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3. 

YBDYMHB 






EL LX 

230.4 

0 . 

445.0 

85.0 

230.4 

85 . 0 

LINE 

445.0 

85.0 

1475.0 

85.0 



CUBI 

1475.0 

85.0 

1822.0 

0.0 

1685 . 0 

85 . 0 

1 . 

YCENTER 






LINE 
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0. 

1822.0 

0. 



10 . 

2BDYTCL 






EL LX 

230.4 

-15.0 

270.0 

40. 0 

230.4 

22 . 0 

LINE 

270 . 0 

40 . 0 

310.0 

70.0 



CUBI 

310.0 

70.0 

400.0 

85.0 

330.0 

85 . 0 

LINE 
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85.0 

650 . 0 

85.0 



CUBI 

650 . 0 

85.0 

785.0 

110.0 

700 . 0 
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CUBI 
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1 050.0 

85. 0 
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LINE 
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85.0 
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LINE 

1250.0 

94 . 0 

1350.0 
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LINE 
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3 . 
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LINE 
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LINE 
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LINE 
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LINE 

794.0 

-85. 0 

900 . 0 

-75.0 



LINE 

900.0 

-75.0 

1100.0 

-75.0 



LINE 
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LINE 
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LINE 
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LINE 
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LINE 
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3. 
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LINE 
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LINE 
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LINE 
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YCENTER 
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YBDYL5CP 

YBDYMHB 






ZBDYLSCP 

ZBDYBCL 





C-141 

INBOARD NACELLE 





0. 

667.91 

883.32 

278.91 

-95.00 

30.0 

0 . 

0,. 

2 . 

0 . 

1 . 

2. 



0 . 

0 .25 

0 .50 

1 . 0 

1 . 5 

2.5 

3.51 

4 . 52 

5.52 

6 . 53 

7.53 

8.54 

9.54 

10.54 

12.55 

14.56 

16.26 

54 . 18 

55.23 

58.75 

63.77 

68.79 

72.81 

77.82 

81.84 

85.86 

89.88 

93.89 

98 . 99 

100.0 






12.53 

12.92 

13.10 

13.37 

13.58 

13.92 

14.18 

1 4 . 40 

14.59 

14.76 

14.90 

15.01 

15.11 

15.18 

15 . 26 

15.31 

15.32 

15 . 32 

15.32 

15.26 

15.06 

14.72 

14.35 

13.75 

13.17 

12.49 

11.72 

10.85 

9.61 

9.34 






C-141 

OUTBOARD NACELLE 





0. 

755.02 

970.43 

456.89 

-95.00 

30. 0 

0 . 

0 . 

1 . 0 

0 . 

1 . 

2. 



0. 

0.25 

0.50 

1 . 0 

1.5 

2.5 

3.51 

4.52 

5.52 

6.53 
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15.26 
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9.34 







85 



OUTPUT DATA ( PRINT £ PLOT) FORMAT 

Printed and plotted output data is provided. Since a typical print or plot 
sequence is lengthy, only a brief description of each type of output will be 
provided here. Note that samples of plotted output can be found in the re- 
sults and geometry verification section of this report. 


Printed Output 

The printed output can be divided into three distinct sections. 

Section I Input Data Listing/Geometry and Grid System 

Verification 

Section II Relaxation Solution Convergence History 

Section III Computed Velocities, Pressures, Forces, Moments, 

Reference Lengths and Areas 

Within each section, the output data will be printed in the following sequence. 


Section I 

• Input Data Listing 

• Case Flow Condition 

• Nominal Extent of Fine Embedded Grid Systems 

• Quick-Geometry Model Error Diagnostic Information (BKOD = ±3 only) 

• Configuration Position in Global Crude Grid System* 

• Body Crude/Fine Grid Limiters (See Note #1, following page) 

• Global Crude Grid Coordinates 

• Fine Embedded Wing Grid Coordinates 

• Pylon Limiters (PY>0) 

• Pod/Nacelle Limiters and Boundary Conditions (Pod> 0) 

• TAG Grid System Coordinates (VER f 0) 

• Winglet Grid System Coordinates (VER>0) 

• Fine Embedded Body Grid Coordinates 


* The configuration is positioned to prevent crude cartesian grid points from 
falling near the wing leading edge. 



Section II 

Phase 1 (Global Crude Grid Solution) 

The following information is printed for each crude grid iteration. 
r • Iteration Count 

• Maximum Correction to the Flow Field Potential (a* ,„ t , r ) 

A / MAX' 

• Grid Position of A$ 

MAX 

Wing Spanwise Circulation (T) Distribution 
NOTE #1 

The body crude and fine grid limiters are the J and K values which 
define special grid points. These limiter values represent the first grid point 
relaxed beyond a fixed potential body boundary point. By using the J and K 
limiters listed, a sketch of the body boundary surface can be composed. A 
typical limiter listing can be found in the table below: 


FINE GRID BODY BOUNDARY POINT LIMITERS 

J 

KUP 

KLO 

JSD 

1 

20 

6 

10 

2 

20 

6 

10 

3 

20 

6 

10 

4 

20 

7 

10 

5 

19 

7 

10 

6 

19 

7 

10 

7 

18 

7 

10 

8 

17 

7 

10 

9 

16 

7 

10 

10 

15 

8 

10 

CRUDE GRID 

BODY BOUNDARY POINT LIMITERS 

J 

KUP 

KLO 

JSD 

1 

22 

15 

4 

2 

21 

15 

4 

3 

20 

15 

4 

4 

19 

16 

4 
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The body boundary surface points have been sketched below in Fig. A- 3. 


CRUDE GRID FINE GRID 



3456789 10 


R84-1 137-035D 


Figure A-3 Body Boundary Point Limiters 


Following the body boundary point limiters, additional limiters will be 
listed if the embedded fine body grid option has been selected. The following 
figure indicates that IBG1, IBGL , JBG , KBGU and KBGL grid points represent 
the first crude point relaxed beyond the embedded fine grid boundary. 



R84-1 13 7-036 D 


Figure A-4 Embedded Body Fine Grid Boundary Limiters 
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The following information is printed at the end of the Global Crude Grid 
Solution . 

r • Section Mach Chart (0 indicates subsonic flow) 

(1 indicates supersonic flow) 

g ) (8 indicates wing section surface) 

(7 indicates wake surface) 

• Span row number/ ti - position/local chord length, chord position, 

V. pressure coefficient, disturbance velocity, section circulation (T). 

Phase 2 (Crude/Fine Grid Interactions) 

The following information is printed for each crude/fine grid iteration 
cycle . 

• Phase 1 (A) Output for Fine Embedded Wing Grid System 

• Pylon Vertical Circulation Distribution (PY>0) 

• Phase I (A) Output For Fine Embedded Winglet Grid System (VER>0) 

• Phase I (A) Output For Fine Embedded Pod Grid System (FIFP > 0 
Only) 

• Phase I (A) Output For TAG (VER i 0) 

• TAG Circulation Distribution For Wing and Winglet 

• Phase 1 (A) Output for Fine Embedded Body Grid System 

(for BKOD > 0 only) 

• Phase 1 (A) Output for Global Crude Grid System 

• Wing Spanwise Circulation (T) Distribution 

For viscous interaction cases, the following information is printed every 
20th cycle. 

• Wing Upper/Lower Surface Boundary Layer Separation Point (x/c) 

• Boundary Layer 6* Slope Added to Wing Boundary Conditions. 

Section III 

The following information is printed at the end of the solution process. 

• Phase 1 (B) Output for Wing Embedded Fine Grid System 

• For (VISMOD = 2,3) Wing Upper/Lower Surface Boundary Layer 
Separation Point (x/c) 
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c c c c 

• Wing Section $- w , nfy, d^, f w 

C 

Note: is section moment about local quarter chord position 

• Span wise Load, Moment, Drag Distributions 

• Wing Exposed Area, Total Area, Aspect Ratio, Taper Ratio, Mean 
Aerodynamic Chord, Average Chord, X-position about which Moments 
are computed 

• Total Wing L WING » M VVING’ D WING 

• Pylon Surface Pressures 

• Phase I (B) Output For Winglet Embedded Fine Grid System 

• Winglet Spanwise Load, Moment and Drag Distributions 

• Winglet Area, Average Chord 

• Winglet C L , C M , C D , d FRICXIO n 

• Pod/Nacelle Pressure Coefficients (see Figure A* 5 for angular 
argument sign conventions) 


90° 



180° 

R84-1 1 3 7-03 7 D 


Figure A-5 Sign Convention for Angular 
Argument of Fuselage/Nacelle 
Pressure Output Stations 
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I 







Pod/Nacelle Surface and Projected Areas 

Pod C L , C M , C D , D FRlCTION 

Pylon Surface Area and Friction Drag 

n 

Wing Pressure Drag D„, . 

P(w) 

Wing Friction Drag C F,., T1 . T _ 

WIN U 

Body Grid Mach Chart 

Body Angular Cut Pressure/Velocity Distribution 
Body Longitudinal Load and Drag Distribution 

Body Length, Wetted Area, Projected Area, Max. Cross-sectional Area, 
Reynolds Number based on body length 
Total Body ^L 


"M 


D, 


BODY’ ‘"BODY’ P(b) 


Body Pressure Drag ^D 


P(b) 


Body Friction Drag r 
Total Configuration C. 


B 


C M’ C D 


Wing-Body Wave Drag C D llIA , rT7 

WAVE 

Wing- Body Friction Drag C„ 

r 

Wing- Body Lift- Induced Drag C Dj 
Wing- Body Spanload Efficiency Factor E 


Plotted Output 

The plot output can be divided into two separate sections. 


Section I (Input Geometry Verification) 

• Title/Case/Flow condition Label 

• Body Cross-Sections 

• Input Wing Sections 

• Configuration Planview 

• Configuration Head-On View 

Section II (Computed Results) 

• Superimposed wing pressure distributions (upper/lower surface) with 
total wing C L> C M , C D label. 

• Wing planform with section shapes at computed span stations. 
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• Detailed wing section pressure distributions with Section C^, C m> C d 
label. 

• Detailed Winglet Pressure Distributions with Section C^, C m> C d label 

• Detailed Pod/Nacelle Pressure Distributions 

• Detailed body angular cut pressure distributions. 

c c c 

• Body longitudinal load plot with body Lg, Mg, Dg and wing-body 
C L> C M > C D label. 

• Wing-body spanload plot with span efficiency, lift-induced drag, wave 
drag, friction drag label. 

INPUT GEOMETRY VERIFICATION 

Coding for graphical inspection of input wing and body geometry has 
been included. Since the program requires considerable computer time and 
core storage to operate and some facility budgets may not provide for a 
number of error- filled submittals, it is recommended that the geometry 
verification mode of operation be used before submitting a complete and 
expensive relaxation solution. The graphic output coupled with printed output 
for geometry and grid systems should be sufficient to diagnose user input 
errors. In particular, it has been found advantageous to make the geometry 
check and perform a single crude and a single fine grid iteration (very cheap) 
to check code flow before a complete analysis is performed. 

The body cross-section array is first to be plotted after case title and 
flow conditions are listed (see Figure A-6). In this case, an error in speci- 
fying a Z-coordinate of the canopy top centerline is apparent. In Figure A-7, 
an error in specifying a Y- coordinate of the windshield base is illustrated. 
Finally, Figure A- 8 illustrates the corrected and final shape of the fuselage to 
be analyzed. 

The defining wing chord sections are displayed after the body geometry 
(see Figure A-9). Each is plotted to a ten inch chord so errors in coordinates 

will become visible. It is important to note that there is no mapping involved 
in the present method; thus, no coordinate smoothing or manipulation is em- 
ployed. Irregularities in input coordinates will cause oscillations in computed 
pressure distributions. 


92 



GULFSTREAM II ISOLATED FUSELAGE 
ISOLATED BODY CASE 1 
MACH = 0.800 ALPHA = 3.10 
RE = 1.00 M 



BODY CROSS-SECTIONS 


R84-1 137-038D 


Figure A-6 Sample Input Geometry Verification Plot with Error in Z-Coordinate 
of Canopy Definition 



Figure A-7 Sample Input Geometry Verification Plot with Error in Y-Coordinate 
of Canopy Definition 
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GULFSTREAM II ISOLATED FUSELAGE 



BODY CROSS-SECTIONS 


R84-1 137-040D 

Figure A-8 Sample Input Geometry Verification Plot for Fuselage 




INPUT WING SECTION 2 2Y/B=0.12 




INPUT WING SECTION 3 2Y/B = 0.35 


R84-1 137-041D 


Figure A-9 Sample Input Geometry Verification Plot for Wing Sections 



A configuration plan -view will follow the wing section figures. This 
should be used to insure that the wing planform is properly defined and its 
placement on the fuselage is correct (see Figure A-10). Finally, a head-on 
view is plotted (see Figure A-ll). Once again, check to see that the wing 
and body are indeed attached. Geometry verification plots which illustrate 
typical nacelle, pylon and winglet positioning checkouts can be seen in Figure 
A-12. 



Fig. A-10 Sample Input Geometry Verification Plot, Plan View 



WING PLANE 


HEAD-ON VIEW 


R84-1 1 3 7-043 D 

Fig. A-11 Sample Input Geometry Verification Plot, Head-On View 
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SUBROUTINE CALL SEQUENCE 


MAIN 



SPLINE 

0 

| — 0 

1 SAREA 


— SLOPP - — 

SPLINE 

oLUrt ■' — 

TAGSL __ 

- orLINfc 

Qp 1 |M c 

■ . SLOPV 

<«« r m | v t 

SPLINE 


NOTE: 


© FUSELAGE GEOMETRY DEFINITION 

© FUSELAQ E geometry interrogation 

RELAXATION/BOUNDARY LAYER ROUTINES 



QWIKDE 


SUBF 


i CSMDEF 


CSMCHK 


BLMDEF 


BLMCHK 


GEMOUT 


CALL SEQUENCE 




DSETUP 


1 DLOKUP 




QWIKLO 


SUBROUTINE CALL SEQUENCE 



CSGEOM 



BLMSET- CURVES 
BLGEOM 


CSMSET 


C BLGEOM 
CSMCOE 
g' "“ ■ " 1 AORDER 


THELIM 


CSMFLT 


CSCALC 


1 — CSMINT 

— SINCOS 

— VDOTV 
MDOTV 


■— SMOTH 


SINCOS 

CSCALC 

BLGEOM 

CSMCOE 

LINLIN 

LINELL 

ELLELL 


SINCOS 

VDOTV 

MDOTV 

CROSS 

ELLCAL 

ELLCAL 

CROSS 

ESTNXT 

SETNXT 
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SUBROUTINE CALL SEQUENCE 



CONTRL- 


■GLOBAL 

•BODVAL 

.PODVAL 

■ INTURP 
•INTERB 

■ POCRUD 

■WINGF - 

■ BODYF 
■BODFIX 
-POFINE 


TRID 

pnn 

tr in 

ruu 

i n i u 

TRID 

TRID 

bodyc “ 

WINGLET 


FILL 


■TRID 

. TAGFIL - TFILL 
-INTURP FILL 

■TRID 

-PLOTER 

-INTEGW 

-INTEG SPLINE 

-INTURP FILL 

- INTURV FILLV 

-LIDRAG MAST - 


0 


PODCP 
SPLINE 
■ SMOTH 
■BODFM- 


PODFM 


-FIT2 

-DELTA1 

- INTRP 

-SERIES - 

- PDAREA 


BLAR 


1 BUDFM — I 

VINTER * — (jO 


— SMOTH 

TAGFIL TFILL 

TAGREL TRID 

INTURV FILLV 

PODF TRID 

VINGF 


c 


INTURV 

TRID 


FILLV 


RFAST 


INTERP 



SUBROUTINE CALL SEQUENCE ( D 


STRIPK 


BLLAM 


VBRAD 


SLOPY 

SPLN1 

SPLN1X 

ORDIN 

VNUSUB 

SLOPBL 

TANCAL 

SPRINT 

FINT 

SOLVEB 

ORDIN 

REDUCX 

RLORD 

GORD 

SIMPSN 


GRAD 

GORD 
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SUBROUTINE DESCRIPTION 


AORDER 

BGLIM 

BLAR 

BLGEOM 

BLLAM 

BLMCHK 

BLMDEF 

BLMSET 

BODFIX 

BODFM 

BODLIM 

BODVAL 

102 
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Orders a set of numbers by permutation index 

Computes wing fixed potential surface limiters in fine embedded 
body grid 

Main control routine for laminar and turbulent modified 
chordwise boundary layer calculation. Computes boundary layer 
displacement thickness (6*) slope for viscous/in viscid interaction 
mode of operation 

Assigns body line model values and derivatives to 
control point coordinates 

Computes Thwaites laminar boundary layer with Rott and Crabtree 
compressibility modification 

Correlates and checks the input data deck and the 
indices for the generated body line math models 

Defines body line models from the input data 

Controls the determination values and first and second derivatives 
for all body line models at a given x-station 

Computes potentials on fixed wing/wake surface in body fine grid 
given solution in fine wing grid and global crude grid 

Computes integrated body force and moment coefficients 

Computes J and K limiters for body boundary in both crude and 
fine grid systems 

Computes body boundary point potential values 



BODYC 

BODYF 

CONTRL 

CROSS 

CSCALC 

CSGEOM 

CSMCHK 

CSMCOE 

CSMDEF 

CSMFLT 

CSMINT 


Finite difference approximations and relaxation solution for body 
boundary in crude grid 

Finite difference approximations and relaxation solution for fine 
body grid 

Main control routine for relaxation solution of governing 
equation, interpolation, boundary layer analysis and 
printed/plotted output 

Solves for the intersection of two lines in a plane 

Computes radial position and derivatives for specified 
cross-section model, arc, and 0' 

Is the main subroutine in the look-up portion of the 
QUICK System. It is called to establish r'= f(0',x). It calls 
appropriate subroutines to evaluate body line values and con- 
struct cross-section geometry at a given x-station. It is used for 
all geometry model interrogation 

Correlates and checks the input data deck and the 
indices for the cross-sectional math model 

Composes the equations which are to define the 
cross-section geometry at a given station 

Logically defines the cross-section models from the input data 

Creates control point definitions to permit the insertion of a 
smooth fillet between cross-sectional arcs 

Locates user specified intersections between cross-sectional arcs 
and adjusts their use- theta limits 
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CSMSET 

CURVES 

DELTA1 

DLOKUP 

DSETUP 

ELLCAL 

ELLELL 

ESTNXT 

FAST 

FILL 

FILLV 

FINT 

FIT2 
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Sets up the control point coordinate arrays used to define the 
cross-section geometry at a specified x -station 

Calculates values and first and second derivatives for 
individual curve fits 

Interpolation routine for wing' spanload 

Is a simple look-up routine. It assigns an index to 

match an input name to a codeword list, but is not capable of 

adding new items to that list 

Is an adapting dictionary look-up routine. New items are added 
to a codeword list, an index (counter) is returned for the 
codeword, and an indicator (INEW) is set equal to 1 when a new 
item is encountered 

Set up for ellipse 

Calculates intersection of two ellipses 

Estimates non-linear root by modified inverse quadratic 

Fast Fourier transform of complex data 
Performs interpolation controlled by INTURP 
Assists INTURV interpolation process 
Simultaneous triple interpolation 

Determines cubic spline fit coefficients for input spanload 
distribution 



GEMOUT 

GEVER 

GLOBAL 

GORD 

GRAD 

INTEG 

INTEGW 

INTERB 

INTERP 

INTRP 

INTURP 

INTURV 

KRVDEF 


Ensures that all body lines required by a cross-sectional 
model are defined for the range of that model 

Controls geometry verification plotting 

Finite difference approximations and relaxation solution 
for global crude grid 

Bradshaw's G function 

Slope of a function at its tabulated points 

Integrates wing load distributions for lift, moment and drag 
coefficients 

Intergrates winglet section coefficients to obtain total lift, moment 
and drag 

Interpolation routine for body fine/global crude grid 
communication 

Interpolation for pod crude/fine grid interaction 

Converts input spanload distribution to a fine over spaced 
distribution 

Controls interpolation for filling fine mesh points using" crude grid 
potential values. Updates crude mesh given fine solution. 

Updates fine mesh given crude solution 

Interpolation for winglet fine/TAG/crude interaction 

Calculates coefficients for the various curve fits 
associated with body line math models 
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LIDRAG 

LINELL 

LINLIN 

MAIN 

MAST 

MDOTV 

ORDIN 

PDAREA 

PLOTER 

POCRUD 

POD 

PODCP 

PODF 

PODFM 

PODSET 


Main control routine for computing' lift induced drag efficiency "E" 
using a Fourier analysis 

Solves for the intersection of a line and an ellipse 
Solves for intersection of two lines 

Controls reading of input data and sets boundary and initial 
conditions. Sets up arrays and storage areas. Sets up crude 
and wing and body fine coordinate systems 

Controls cubic spline fit for interpolation of input 
spanload distribution 

Performs matrix multiplication of a vector 
Linear interpolation 

Computes pod surface and projected areas 

Controls all graphic output (except input geometry verification) 

Prints results in global crude grid; mach chart, wing 
upper/lower surface pressure coefficients 

Pod crude grid relaxation process 

Compute pod/nacelle pressure coefficients 

Pod/nacelle fine grid relaxation process 

Computes pod/nacelle forces and moments 

Set up for pod/nacelle simulation. Pod geometry input 



PODVAL 

POFINE 

PYSET 

QWIKDE 

QWIKLO 

REDUCX 

RFAST 

RLORD 

SAREA 

SERIES 

SETNXT 

SIMPSN 

SINCOS 

SLOPP 

SLOPV 

SLOPBL 


Computes pod/nacelle boundary condition for crude and 
fine grid 

Prints output results for fine grid arrays 

Set up for pylon surface simulation. Pylon geometry input 

Main control routine for Quick- Geometry definition and check out 

Main control routine for interrogation of Quick- Geometry math 
model 

Performs interpolation to new grid 
Fast Fourier transforms of real data 
Bradshaw's L function 

Computes body surface area given an array of cross-sections 
Determines Fourier series coefficients 
Reorders points for non-linear root finder 
Simpson's rule integration 

Adjusts input interrogation angles for top and bottom dead center 

Computes pylon boundary conditions 

Computes winglet boundary conditions 

Slope of a tabulated function at an arbitrary point 
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SLOPE 

SLOPY 

SMOTH 

SOLVEB 

SPLINE 

SPLN1 

SPLN1X 

SPRINT 

STRIPK 

TAGFIL 
TAG SET 
TAGSL 
TAGREL 
TANCAL 

TFILL 


Computes boundary conditions for wing- surface and axisymmetric 
bodies 

Computes wing surface slopes 
Function for smoothing an array of values 
Solution of two simultaneous linear algebraic equations 
Computes a cubic spline through a set of points 

Computes continuous derivatives interpolation by means of a cubic 
fit 

Entry for special cases requiring extrapolation beyond ends of X 
and Y tables 

Prints output of profile results 

Starting condition setup and flow control for laminar/turbulent 
boundary layer prediction 


Interpolation for TAG/crude grid interaction 


Set up for wing tip augmentation grid TAG 
Computes wing boundary conditions for TAG 
Relaxation process for TAG 

Computes characteristic angles for use in B.L. solution 
(equation 21 of Bradshaw and Ferriss) 

Assists TAGFIL interpolation process 
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THELIM 

TRID 

VBRAD 

VDOTV 

VERSET 

VINGF 

VINTER 

VNUSUB 

WINGF 

WINGLET 


Creates and controls use- theta arrays to establish continuity in 
the cross-section model 

Solves tri-diagonal matrix 

Computes Bradshaw compressible 2-D turbulent boundary layer 
simulating' 3-D boundary layer on infinite yawed wing by 
Nash-Tseng modified chord technique 

Computes a vector dot product 

Set up for winglet fine embedded grid system. Winglet geometry 
input 

Winglet fine embedded grid relaxation process 

Performs cubic 6* fit for separated boundary layer in wing 
section cove regions 

Computes the Nash effective viscosity 

Finite difference approximations and relaxation solution for wing 
fine grid 

Crude grid relaxation process for winglet in crude grid 
system 


109 



KEY VARIABLE DESCRIPTIONS 

This description of key program variables and constants which are located 
in several common blocks will be useful in understanding' flow logic. 

VARIABLE D ESCRIP TION 

AK The value l-M^ 

ALAM Wing taper ratio (X) 

ALPHA Angle-of-attack (radians) 

AMAC Wing mean aerodynamic chord (MAC) 

AMACH Mach number 

AM2 The value M^ 

AOA Angle of-attack (degrees) 

AR Wing aspect ratio (AR) 

BAREA Body wetted area 

BCF Body skin friction coefficient 

BCL Wing crude grid lower boundary slopes 

BCLF Wing fine grid lower boundary slopes 

BCLP Pylon outboard boundary condition 

BCLT Wing lower surface boundary condition in TAG 

BCLV Winglet outboard boundary condition 
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VARIABLE 


DESCRIPTION 


BCU 

BCUF 

BCUP 

BCUT 

BCUV 

BNOSE 

BNOSEP 

BODCD 

BODCL 

BODCM 

BP ARE A 

BS 

BTAIL 

BTAILP 

CA 

CAV 


Wing crude grid upper boundary slopes 

Wing fine grid upper boundary slopes 

Pylon inboard boundary condition 

Wing upper surface boundary condition in TAG 

Winglet inboard boundary condition 

X coordinate of body nose 

X coordinate of pod/nacelle nose (input as BNOSE 1) 

Body (integrated) drag coefficient 

Body (integrated) lift coefficient 
Body (integrated) moment coefficient 

Body projected area 

Body plot scaling coefficient 

X coordinate of body tail 

X coordinate of pod/nacelle tail (input as BTAIL 1) 
Global crude grid stretching coefficient (£ ) 

X 

Wing average chord (C A . ) 

A V 

Winglet average chord 


CAVW 


VARIABLE 

DESCRIPTION 

CB 

Global crude grid stretching' coefficient (£ ) 

cc 

Global crude grid stretching coefficient (n ) 

CCT 

TAG stretching coefficient (ti ) 

CD 

Global crude grid stretching coefficient (n ) 

CDI 

Lift induced drag coefficient 

CDINT 

Wing section integrated drag 

CDINW 

Winglet section drag coefficients 

CDT 

TAG stretching coefficient (ri ) 

CE 

Global crude grid stretching coefficient (c ) 

CET 

TAG stretching coefficient (t ) 

CF 

Global crude grid stretching coefficient (c ) 

Z Z 

CFINT 

Integrated wing section friction coefficient 

CFT 

TAG stretching coefficient (c ) 

CIR 

Wing circulation (T) 

CIRPC 

Pylon circulation at crude grid vertical cuts 

CIRPY 

Pylon circulation at fine grid vertical cuts 

CIRV 

Winglet circulation at vertical span cuts 



VARIABLE 

CIRVC 

CLINT 

CLINW 

CMINT 

CMINW 

CMLOC 

CMLOW 

CPL 

CPU 

CPUPOD 

CSCUT 

DELSL 

DELSU 

DETA 

DIM 

DRDXC 


DESCRIPTION 

Winglet circulation at crude grid vertical cuts 
Wing section integrated lift 
Winglet section lift coefficients 
Wing section integrated moment 
Winglet section moment coefficients 

Wing section integrated moment about local quarter chord 

Winglet section local moment coefficients 
Wing lower surface pressure coefficient 

Wing upper surface pressure coefficient 

Pod pressure coefficients 

Body x station for cross-sectional cut 

Wing boundary layer slopes for section lower surface 

Wing boundary layer slopes for section upper surface 

Global crude grid mesh spacing in n direction 

Configuration length for non-dimensionalizing maximum potential 
updates 

Axisymmetric body slope distribution in crude grid 
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VARIABLE 


DESCRIPTION 


DRDXCP 

DRDXF 

DRDXFP 

DXB 

DXBP 

DXI 

DXV 

DXW 

DYB 

DYBP 

DYV 

DYW 

DZB 

DZBP 

DZETA 

DZW 


Pod/nacelle boundary condition 

Axisymmetric body slope distribution in fine grid 
Pod fine grid boundary condition 
Fine body grid mesh spacing in X direction 
Fine pod grid mesh spacing in X direction 
Global crude grid mesh spacing in £ direction 
Winglet fine grid mesh spacing in X direction 
Fine wing grid mesh spacing in X direction 
Fine body grid mesh spacing in Y direction 
Fine pod grid mesh spacing in Y direction 
Winglet fine grid mesh spacing in Y direction 
Fine wing grid mesh spacing in Y direction 
Fine body grid mesh spacing in Z direction 
Fine pod grid mesh spacing in Z direction 
Global crude grid mesh spacing in c direction 
Fine wing grid mesh spacing in Z direction 



VARIABLE 


DESCRIPTION 


E 

ETA 

ETAT 

FMFR 

FNPR 

G 

H 

IBGI 

IBGIP 

IBGL 

IBGLP 

IL 

ILEF 

ILEV 

ILT 

ILV 

ILVC 


Wing spanload efficiency 

ti coordinates for global crude grid (transformed space) 

TAG 11 coordinates 
Nacelle inlet mass flow ratio 
Nacelle exhaust nozzle pressure ratio 
The value (JT+ 1 )M^ 

The value (y-l)M^ 

Crude grid I value of body grid inner overlap region forward boundary 
Crude grid I value at pod fine grid forward boundary 
Crude grid I value of body grid inner overlap region aft boundary 
Crude grid I value at pod fine grid aft boundary 
Crude grid wing leading edge I value 
Wing fine grid leading edge I value 
Fine grid winglet leading edge I value 
TAG wing leading edge I value 
TAG winglet leading edge I value 
Crude grid winglet leading edge I value 
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VARIABLE 


DESCRIPTION 


IMACH 

IMAX 

IMAXB 

IMAXBP 

IMAXV 

IMAXW 

INOSE 

INOSEC 

INOSEP 

INOSPF 

IPTYPE 

IT 

ITAGI 

ITAGIT 

ITAGL 

ITAIL 

ITAILC 


Code for subsonic (0) or supersonic (1) flow at a grid point 

Maximum number of crude grid points in X direction 

Maximum number of fine body grid points in X direction 

Maximum number of fine pod grid points in X direction 

Maximum number of winglet fine grid points in X direction 

Maximum number of fine wing grid points in X direction 
Fine body grid I value at body nose 

Crude grid I value at body nose 

Crude grid I value at pod nose 

Fine grid I value at pod nose 

Nacelle type - closed body/hot jet/cold jet 

Crude grid wing trailing edge I values 

Crude grid I limiter for forward TAG boundary 

On/off code for wing tip augmentation grid (TAG) 

Crude grid I limiter for aft TAG boundary 

Fine body grid I value at body tail 

Crude grid I value at body tail 



VARIABLE 


DESCRIPTION 


ITAILP 

ITALPF 

ITEF 

ITER 

ITEV 

ITMAX 

ITT 

ITV 

ITVC 

JBG 

JBGIP 

JBGLP 

JBY 

JMAX 

JMAXB 

JMAXBP 

JMAXV 


Crude grid I value at pod tail 
Fine grid I value at pod tail 
Wing fine grid trailing edge I value 
Iteration count 

Fine grid winglet trailing edge I value 
Maximum number of TAG grid points in X direction 

TAG I value of wing trailing edge 

TAG I value of winglet trailing edge 

Crude grid I value of winglet trailing edge 

Crude grid J value of body grid inner overlap region side boundary 

Crude grid J value at pod fine grid inboard boundary 

Crude grid J value at pod fine grid outboard boundary 

Fine grid J value of winglet section plane 

Maximum number of crude grid points in Y direction 

Maximum number of fine body grid points in Y direction 

Maximum number of fine pod grid points in Y direction 

Maximum number of winglet fine grid Y-stations 
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VARIABLE 


DESCRIPTION 


JNOSEP 

JROOT 

JSD 

JSDC 

JSDPL 

JSDPR 

JTAGI 

JTAGO 

JTIP 

JTIPT 

JTMAX 

KBB 

KBC 

KBGL 

KBGLP 

KBGU 

KBGUP 


Crude grid J value at pod position 
Crude grid J value at wing root 

Fine body grid J value at first influence of body boundary points 

Crude grid J value first influenced by body boundary 
point 

Fine pod grid J limiter at pod side-left 
Fine pod grid J limiter at pod side-right 
Crude grid J limiter for inboard TAG boundary 
Crude grid J limiter for outboard TAG boundary 
Crude grid J value at wing tip 
TAG J value at wing tip 

Maximum number of TAG grid points in Y direction 
Fine body grid K value at wing plane 
Crude grid K value at wing plane 

Crude grid K value of body grid inner overlap region (lower) 
Crude grid K value at pod fine grid lower boundary 
Crude grid K value of body grid inner overlap region (upper) 
Crude grid K value at pod fine grid upper boundary 



VARIABLE 


DESCRIPTION 


KBOD 

KBW 

KLOC 

KLOF 

KLOPF 

KMAX 

KMAXB 

KMAXBP 

KMAXW 

KNOSEP 

KODB 

KOPOD 

KOPY 

KOVER 

KPLO 

KT 


Code for body in crude grid (0) or fine embedded grid (1) 

Fine wing grid K value at wing plane 

Crude grid K limiters for body surface (lower) 

Fine grid K limiters for body surface (lower) 

Nacelle fine grid lower surface K limiters 

Maximum number of crude grid points in Z direction 

Maximum number of fine body grid points in Z direction 

Maximum number of fine pod grid points in Z direction 

Maximum number of fine wing grid points in Z direction 

Crude grid K value at pod position 

Body option code. .. cylinder, axisymmetric, arbitrary 
body (input as BKOD) 

Number of pods/nacelles (input as PCD) 

Number of pylon surfaces (input as PY) 

Number of winglets (input as VER) 

Crude grid lower K limit of pylon surface 

TAG K value of wing plane 


VARIABLE 


DESCRIPTION 


KTAGL 

KTAGU 

KTC 

KTIPV 

KTMAX 

KUPC 

KUPF 

KUPPF 

MAXIT 

MAXITF 

MODV 

NCASE 

NINB 

NINBP 

NINP 

NINV 


Crude grid K limiter for lower TAG boundary 

Crude grid K limiter for upper TAG boundary 

Crude grid K value at winglet tip 
TAG winglet tip K value 

Maximum number of vertical (Z) points in TAG 
Crude grid K limiters for body surface (upper) 

Fine grid K limiters for body surface (upper) 

Nacelle fine grid upper surface K limiters 

Maximum number of initial crude grid iterations (input as AXIT) 

Maximum number of crude/fine grid cycles (input as AXITF) 

Mode of operation for viscous effects (input as VISMOD) 

Case description. . .wing, body or wing- body case (input as 
CASE) 

Number of ordinates defining axisymmetric body shape 
(input as BNIN) 

Number of pod ordinates defined (input as BNIN 1) 

Number of ordinates defining each pylon section (input as PIN) 
Number of ordinates defining each input winglet section 


(input as VIN) 



VARIABLE 


DESCRIPTION 


NINW 

NOSEB 

NOSEBP 

NOSEW 

NOSP 

NOSW 

NPOA 

NPOB 

NPSEC 

NSECT 

NTC 

NTF 

NTOTAL 

NTOTB 

NTOTV 


Number of ordinates defining each wing section (input as ANIN) 

Blunt/sharp nose body code (for spline fit) 

Sharp/blunt nose shape code for pods 

Blunt/sharp nose wing code (for spline fit) 

Pylon sharp/blunt nose code 

Sharp/blunt winglet section nose shape code 

Number of fine X-grid points between leading and trailing 
edge of each wing section 

Number of fine body X-grid points between nose and tail of 
body 

Number of defining pylon sections (input as PSEC) 

Number of defining wing sections (input as ASECT) 

Number of points representing body cross-sections in crude grid 

Number of points representing body cross-sections in fine body 
grid 

Total number of crude grid points in single X-Z plane 
Total number of fine body grid points in single X-Z plane 
Total number of fine winglet grid points in single X-Y plane 
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VARIABLE 


DESCRIPTION 


NTOTW 

NVSEC 

NWPO 

PAREA 

PBL 

PCI 

PC2 

PC3 

PCF 

PCL 

PF1 

PF2 

PF3 


Total number of fine wing grid points in single X-Z plane 

Number of winglet sections defined (input as VSEC) 

Code for print out of crude grid results for diagnostic pur- 
poses 

Pod surface area 

Fine body grid wing/wake lower surface potentials 

Global crude grid potential (<^) arrays 

(Note: Only three planes are in core at one time) 

Global crude grid potential (tp) arrays 

(Note: Only three planes are in core at one time) 

Global crude grid potential (</>) arrays 

(Note: Only three planes are in core at one time) 

Pod/nacelle skin friction coefficient 

Crude grid wing/wake lower surface potentials 

Fine embedded wing and body potential arrays. 

(Note: Only three planes are in core at one time) 

Fine embedded wing and body potential arrays. 

(Note: Only three planes are in core at one time) 

Fine embedded wing and body potential arrays. 

(Note: Only three planes are in core at one time) 



VARIABLE 


DESCRIPTION 


PFL 

PI 

PLTAG 

PODALF 

PODBET 

PODCD 

PODCL 

PODCM 

PPAREA 

PTAG 

PVCL 

PVL 

PVTL 

PYAREA 

PYLO 

PYLOC 

RADIUS 


Fine wing grid lower surface potentials 
n 

TAG wing lower surface potential array 
Pod/nacelle incidence relative to fuselage 
Pod/nacelle yaw angle 
Pod/nacelle drag coefficient 
Pod/nacelle lift coefficient 
Pod/nacelle moment coefficient 
Pod projected area 
TAG field potential (0) array 

Winglet outboard surface potential (0) array for crude grid 
Winglet fine embedded grid outboard surface potential (0) array 
Winglet TAG outboard surface potential (0) array 
Pylon surface areas 

Pylon fine grid outboard surface potential (0) array 
Pylon crude grid outboard surface potential (0) array 
Radius of body cylinder (input option) 
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VARIABLE 


DESCRIPTION 


RAVC Crude grid average body radius for boundary condition 

calculation 

RAVF Fine grid average body radius for boundary condition 

calculation 

RC Axisymmetric body radius distribution in crude grid 

RCP Radii for pod surface after positioning 

RE Freestream Reynolds number 

RF Axisymmetric body radius distribution in fine grid 

RFP Pod radii in fine grid 

RIN Radial coordinates defining axisyirmetric body (input option) 

RINP Radial coordinates defining axisyrrmetric pod 

RMAX Body maximum radius for computational body surface 

RMAXP Pod maximum radius 


Wing spanwise drag coefficient CC^ 

~c77 


SDW Winglet spanwise drag coefficient CC d 



SEXP Wing exposed area 
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VARIABLE 

SFD 

SGRAD 

SLD 

SLW 

SMD 

SMW 

TCIR 

TDX 

TDY 

TDZ 

THETC 

THETF 

THETFP 

TITLE 

TLE 

TPIST 


DESCRIPTION 


Wing spanwise friction coefficient CC^. 


C 


AV 


Body side slope at wing- body juncture 
Wing span load coefficient CC^ 


Winglet spanwise load coefficient CC„ 


'AV 


Wing spanwise moment coefficient CC 


m 


'AV 


Winglet spanwise moment coefficient CC 

m 
C AV 

Wing circulation at TAG wing span cuts 

TAG A X spacing 

TAG A Y spacing 

TAG A Z spacing 

Body crude grid angular cuts 

Body fine grid angular cuts 

Angular cuts for pod boundary points 

Case title for identifying graphic and printed output 

Wing leading edge X value at TAG span cuts 

Pylon yaw incidence or twist 
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VARIABLE 


DESCRIPTION 


TSLOC 

TTE 

TVIST 

TWIST 

VANGL 

VAREA 

VCD 

VCL 

VCLE 

VCM 

VCTE 

VLE 

VTE 

W 

WAREA 

WCD 
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Wing local sweep angle at wing fine grid boundary points 

Wing trailing edge X value at TAG span cuts 

Winglet section incidence angle or twist 

Wing twist (incidence) distribution 

Winglet cant angle, measured from wing plane 

Winglet area 

Winglet drag coefficient 

Winglet lift coefficient 

Crude grid X value of winglet leading edge 

Winglet moment coefficient 

Crude grid X value of winglet trailing edge 

Winglet leading edge X value after positioning 

Winglet trailing edge X value after positioning 

Relaxation factor w 

Wing area S^ 

Wing drag coefficient 



VARIABLE 

WCF 

WCL 

WCM 

WCORD 

WS 

X 

XBF 

XBFP 

XI 

XILE 

XINB 

XINEP 

XINP 

XINV 

XINW 

XITE 

XLE 

XLET 


DESCRIPTION 
Wing friction drag coefficient 

Wing lift coefficient 
Wing moment coefficient 
Wing section local chord length 
Wing plot scaling coefficient 

X coordinate for global crude grid (physical space) 

X coordinate for body fine grid 
Fine pod grid X array 

5 coordinate for global crude grid (transformed space) 
i coordinate of local wing section leading edge 
X coordinates defining axisynmetric body (input option) 

X coordinates defining axisynmetric pod 

X ordinates defining pylon 
X ordinates defining winglet section 
X ordinates defining wing section 
£ coordinate of local wing section trailing edge 
X coordinate of local wing section leading edge 
X coordinate of wing tip leading edge 


VARIABLE 


DESCRIPTION 


XMOM 

XNC 

XNCP 

XNF 

XNOSE 

XOL 

XOLP 

XPL 

XPT 

XSF 

XTAG 

XTAIL 

XTE 

XTET 

XVF 

XVL 
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Position about which configuration moments are computed 

Body normal vector X direction at crude grid body points 

Pod normal vector X direction at crude grid points 

Body normal vector X direction at fine grid body points 

X coordinate of body nose repositioned in crude grid 

Non-dimensional distance along body length and wing chord 

X stations for pod surface after positioning 

X coordinate of input wing section leading edge 

X coordinate of input wing section trailing edge 

X coordinate defining wing section at each fine grid plane 
points 

TAG grid X array 

X coordinate of body tail repositioned in crude grid 
X coordinate of local wing section trailing edge 
X coordinate of wing tip trailing edge 
Winglet fine grid X array 

X value of winglet leading edge at input stations 



VARIABLE 


DESCRIPTION 


XVT 

XWF 

Y 

YBF 

YBFP 

YBODYP 

YINL 

YINLP 

YINLV 

YINU 

YINUP 

YINUV 

YNC 

YNCP 

YNF 

YOB 

YP 


X value of winglet trailing edge at input stations 
Fine embedded wing grid X coordinate 

Y coordinate for global crude grid (physical space) 

Y coordinate for body fine grid 
Fine pod grid Y array 

Y coordinate of pod/nacelle axis (input as YBQDY 1) 

Y ordinates defining lower wing section 

Y ordinates defining pylon outboard surfaces 

Y ordinates defining winglet outboard surface 

Y ordinates defining upper wing section 

Y ordinates defining pylon inboard surfaces 

Y ordinates defining winglet inboard surface 

Body normal vector Y direction at crude grid body points 
Pod normal vector Y direction at crude grid points 
Body normal vector Y direction at fine grid body points 
Wing span station (2Y/b) 

Y coordinate of input wing section trailing edge 



VARIABLE 


DESCRIPTION 


YSF 

YTAG 

YTIP 

YV 

YVF 

Z 

ZBF 

ZBFP 

ZBODYP 

ZETA 

ZETAT 

ZNC 

ZNCP 

ZNF 

ZOB 

ZTAG 

ZWF 

ZWING 


Y coordinate defining wing section at each fine grid plane 
TAG grid Y array 

Y coordinate of wing - tip 

Y value of winglet input stations 
Winglet fine grid Y array 

Z coordinate for global crude grid (physical space) 

Z coordinates for body fine grid 
Fine pod grid Z array 

Z coordinate of pod/nacelle axis (input as ZBGDY1) 

C coordinate for global crude grid (transformed space) 

TAG c coordinates 

Body normal vector Z direction at crude grid body points 
Pod normal vector Z direction at crude grid points 

Body normal vector Z direction at fine grid body points 

Winglet span station (y/b) 

TAG grid Z array 

Fine embedded wing grid Z coordinate 
Wing height relative to center of body 
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APPENDIX B 


USAGE/MODIFICATIONS FOR LaRC CDC/VPS-32 VERSION 


The LaRC VPS-32 version of the WIBCO-PPW code, using VSOS 2.1.5 (opt = E) 

requires less than 1 million^g storage for load and execution. Differences 
are noted as follow: 

1. FILES USAGE 


Thirteen (13) disk units are used, including one for a solution 
save/restart. The planar potential arrays are set up and accessed as 
table lookup in a large labeled cannon block, eliminating disk units 2, 

3, 4, and 5 (units 12, 13, 14, and 15 of IBM code) by duimy READMS, 

WRIIMS referencing. Several options are available for the solution 
save/restart capability, using unit 7 (unformatted data). Data needed by 
a separate plot program are saved on units 80, 82, 83, 87, 88 (all 
unformatted) and unit 9 (formatted). 

2. UNIT 1 DATA INPUT 


Through use of unit 1 data input the user can redefine certain key 
parameters that have been code-defined (hard coded). These code-defined 
values should be used initially to enhance the probability of proper 
execution for input geometry. Parameter values revert to default 
settings if keys for changes are not activated. Unit 1 data referenced 
in the information which follows use 8-field 10-digit forrat (8F10.0) 
with decimal point required. Default values for variables are included 
in the description for variables. The changes allowed are: 

a) Number of fine x-grid points at wing tip is specified on Chrd 1-W 
rather than Card 2-W. 

b) Save/restart of solution in the crude/fine grid only. 

c) Redefinition of wing fine grid boundaries. 

d) Redefinition of transition location at all span stations for viscous 
flow. 

e) Redefinition of relaxation factors and interaction frequency in the 
general solution process. 

The notation of appendix A for configuration indentification is used to 
indicate the card image modifications and additions for these allowed 
changes. Note that only the variables redefined from appendix A are included 
here. 
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CARD 

NUMBER 


CARD VARIABLE 
COLUMN NAME 


DESCRIPTION 


Qird 2-A 1-10 CASE |CASE| Key describing the type of 

configuration. 

CASE = 1. Isolated Body (omit cards -W). 

CASE = 2. Isolated Wing (omit cards -B). 

CASE = 3. Wing-Body. 

SGN (CASE) Trigger to control relaxation 
factors and interaction 
frequencies in general solution 
process . 

CASE > 0. Same as original input, no 

additional input required on 
Card 2-A-2 . 

CASE < 0. Read additional input, 

Card 2-A-2. 


61-70 VISMOD I VISMOD I Key describing extent of viscous 


ef f ects . 

VISMOD = 1. 

No viscous effects. 

VISMOD = 2. 

Viscous effects computed at 
end of inviscid analysis. 

VISMOD = 3. 

Inviscid/viscous inter- 
action. 


SGN (VISMOD) Trigger to redefine transition 
location, x/c, for viscous 
flow. 

VISMOD > 0. Same as original code — 
no additional input 
required on Card 6-W. 

VISMOD < 0. Read additional input, 

Card 6-W. 
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CARD VARIABLE 

COLUM N NAME DESCRIPTION 

71-80 START START = 0. Conventional independent run 

(default). 

START = 1. Save final solution at end 
of current AXITF iteration 
cycles. 

START = 2. Restart solution at previous 

AXITF+1 iterations; do not save 
final solution. 

START = 3. Restart solution at previous 
AXITF+1 iterations; save final 
solution at end of current AXITF 
iteration cycles. 

START = 4. Save solution at end of a 

specified crude/fine iteration 
cycle during a conventional 
independent run. This option 
requires additional data, 

Card 2-A-l. 

NOTE: Card 2-A-l required only for START = 4. 

Card 2-A-l 1-10 AXITF2 Crude/fine iteration cycle at which 

solution is saved during a conventional run. 
(AXIT+1 < AXITF2 < AXITF) 

NOTE: Card 2-A-2 required only for CASE < 0. 

Card 2-A-2 1-10 RFSUB Subsonic relaxation factor (default = 1.5). 

11-20 RFBL Relaxation factor for boundary layer 

(default = 0.6). 

21-30 XITBLI Crude/ fine iteration cycle at which boundary 

layer is initially relaxed (default = 5.0). 

31-40 XITBLF Frequency of boundary-layer relaxation in 

crude/fine solution (default = 20.0). 

41-50 XITFBB Frequency of fine body (pod) boundary value 

potential update in crude/fine solution 
(default = 10.0). 

(Note: Any variable set = 0.0 is assigned default value.) 


CARD 

NUMBER 

Card 2-A 
(contd) 
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CARD CARD VARIABLE 

NUMBER COLUMN NAME DESCRIPTION 


C&rd 1-W 1-10 ASECT |ASECT| Number of streamwise sections 

defining wing planform 
(1 < ASECT < 20). 

SGN (ASECT) Trigger to redefine wing fine 
grid boundaries 

ASECT > 0. Same as original code, no 
additional input required 
on Card 1-W-l. 

ASECT < 0. Read additional data, 

Card 1-W-l. 


71-80 AQXFG AQXFG = -1. Code sets 100 fine x-grid points 

at wing tip station (default). 

AQXFG = 0. Code determines the number of 
fine x-grid points at wing tip 
station. 

AQXFG = other. Number of fine x-grid 
points at wing tip 
station. 

(10. < AQXFG < 99.) 


NOTE: 

Card 1-W- 

-1 required only for ASECT < 0. 

Card 1-W-l 

1-10 

AWLE 

Percent of wing fine grid in front of leading 
edge (default = 20.). 


11-20 

AWTE 

Percent of wing fine grid behind trailing 
edge (default = 10.). 


21-30 

AWLS 

Percent of wing fine grid below wing 
(default = 10.). 


31-40 

AWUS 

Percent of wing fine grid above wing 
(default = 30.). 

(Note: 

Any variable set 

= 0.0 is assigned default value. ) 
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CARD 

NUMBER 

Card 2-W 


NOTE: 
Card 6-W 


CARD VARIABLE 

COLUMN NAME DESCRIPTION 


00 Field unused. 

Card 6-W required only for VISMOD < 0. 

1-80 ATRIP Transition location, x/c, at every wing 

span station for viscous flow. Specify 
from inboard to outboard for 18 upper 
surface values, 18 lower surface values, 
respectively. (Default values assigned, 
.05, at all wing span stations.) 


135 


136 



APPENDIX C 

GENERAL RECOMMENDATIONS FOR CODE USAGE 


The computer program described herein has been implemented by a num- 
ber of investigators. To date, over 75 different configurations (see Table 
C-l) have been modeled and analyzed. Because of problem three-dimension- 
ality and inherent code flexibility, however, it is impossible to verify the 
computer program for all of the possible option combinations, shape variations 
or flow conditions that might be of interest. Despite this, good results should 
be obtained if the user is aware of the following recommendations. 

1. PRACTICE 

There are several factors that combine to foil initial attempts by new 
users performing transonic aircraft analysis. These factors are listed below. 

• Complex geometry input/definition 

• Complexity/sensitivity of transonic flow fields 

• Character of finite -difference relaxation process. 

Experience indicates that input and checkout of a complex geometric shape 
such as a fuselage can be a difficult task when setting up for a subsonic panel 
code analysis. For transonic analyses, the problem is aggravated further by 
the complexity of the flow field and sensitivity to what may appear as minor 
surface irregularities. Experience with subsonic "panel" methodology does not 
prepare the user for transonic code applications. 

The problem is compounded by the nature of finite- difference relaxation 
processes which require a computational grid system. A surface definition 
discrepancy may, in addition to causing erroneous computed pressures, result 
in generation of an improper grid system. 

As a result, it is recommended that the user become familiar with the 
computer code by first analyzing simple shapes such as isolated wings and 
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AIRFOILS (2— D) 

NACA 0012 
NASA LS(1 ) 

NASA TRANSPORT AIRFOIL 

WINGS 
ONERA M-6 
RAE 101 

NACA WING (RM A9K01) 

LOCKHEED/AFOSR (3 CONFIGS) 

LTV SKEWED CAMBER WING 

ISOLATED BODIES 

NACA SHARP NOSE BODY (RM L53H04, 2 CONFIGS) 

NASA BLUNT NOSE BODY (TN D-7331 ) 

GRUMMAN GULFSTREAM II FUSELAGE 

STORE CALIBRATION BODY 

DOUGLAS A -4 FUSELAGE 

NIELSEN STORE 

E A-6B TAIL POD 

NTF 5° CALIBRATION CONE 

PATHFINDER I FUSELAGE 

NACA BODY OF REVOLUTION (RM A9I26) 

WING-BODY/WING-FUSELAGE 
GRUMMAN F-1 4 (3 CONFIGS} 

NACA RESEARCH MODEL (RM L51F07) 

NASA RESEARCH MODEL (TM X-3431) 

GENERAL DYNAMICS F-16 
GENERAL DYNAMICS F-1 1 1/TACT 
LAVI FIGHTER 

LTV A-7 SUPERCRITICAL WING CONFIG 
DOUGLAS TRANSPORT MODEL 

GRUMMAN X-29A FORWARD SWEPT WING DEMONSTRATOR* 

BOEING 747-200 
BOEING 747-SP 

NASA F-8 RESEARCH CONFIG 
ROCKWELL HiMAT CANARD DESIGN 
NASA USB TRANSPORT 

GRUMMAN/AIR FORCE STAC/CD AF CONFIGS 
BOEING B-52 

NASA TAILOR-MATE (2 CONFIGS) 

NASA ENERGY EFFICIENT TRANSPORT 

ROCKWELL CANARDED RESEARCH FIGHTER CONFIG. 

GRUMMAN ATF-RESEARCH FIGHTER CONFIG 
BEECH TRANSPORT 

NACA SKEWED WING RESEARCH MODEL (RM A58C03) 

NTF-PATHF INDER II FIGHTER 

COMPLEX CONFIGURATIONS (NACELLES/PODS, PY LONS, WING LETS} 
SPACE SHUTTLE (LAUNCH CONFIG) 

LOCKHEED L-1011 
GULFSTREAM AEROSPACE G-lll 
BOEING KC-135 

DOUGLAS WING-BODY-NACELLE-PROP SLIPSTREAM 

LOCKHEED C-141 

BOEING 767 WINGLET STUDY 

DOUGLAS DC-10 WINGLET STUDY (2 CONFIGS} 

LOCKHEED C-5A 

BOEING 747 WINGLET STUDY 

GRUMMAN A-6F PYLON STUDY 

GRUMMAN DESIGN 698 TILT-NACELLE V/STOL 

FUJI TRAINER (3 CONFIGS} 

GRUMMAN VTX WING/WINGLET DESIGN 
GRUMMAN E-2 WINGLET STUDY 
ARA MODEL M-64 RESEARCH CONFIGURATION 
NASA EET NACELLE/PYLON STUDY (2 CONFIGS) 

ROCKWELL HIMAT 

NTF-PATHF INDER I TRANSPORT 

GULFSTREAM AEROSPACE G-IV WING 

GRUMMAN DESIGN 623 V/STOL AERO-PROPULSION INTEG STUDY 
NASA POD-IN-WING CONFIG 
NASA PROP-FAN SIMULATOR 

R84-1 137-046D 


Table C-1 Configurations Analyzed Using Transonic Wing-Body Code 



isolated bodies at transonic conditions. Complex aircraft shapes should be 
attacked after a basic knowledge of 1) the computer code, and 2) the charac- 
ter of transonic flow fields is in hand. Then, when a problem surfaces, the 
user can determine whether the computer code is at fault, the input data is in 
error, or the flow condition is unreasonable. 

2. NON-POTENTIAL PHENOME NA 

In general, most aircraft operate in a low-gradient or "small disturbance" 
world. Flight efficiency degrades rapidly if disturbances are large. As a 
result, the present small- disturbance formulation has been found to be applica- 
ble to many different configurations over an extraordinary range of flow 
conditions. The user should be aware, however, that there are many cases 
for which the code can not be expected to provide a good flow simulation. If 
geometry and flow conditions combine to produce severe flow gradients, flow 
separation will occur. If the separation region is large (i.e., ~10% chord on 
the wing) the computational and physical flows will not agree. In addition, if 
vortices generated at the wing leading edge or by configuration components 
interact with surfaces that are of interest, flow predictions will be com- 
promised . 

3. DRAG PREDICTION 

Consistent prediction of absolute and incremental drag values is not 
possible with today's 3-D computational methods. This program is not the 
exception. Some useful incremental drag predictions have been obtained 
during project applications, but good results have not been obtained with 
sufficient regularity to recommend usage. In general, the simpler the shape, 
the higher the probability that drag predictions will be useful. 

4. LIFT CONVERGENCE 

All of the 3-D computations included in this report were generated using 100 
crude grid iterations and 80 cycles in the fine/crude system. Experience 
indicates that the resulting convergence level is sufficient for most engineering 
applications. In other words, the accuracy vs cost ratio is high. For some 
high aspect ratio transport wings, however, this level of convergence may 
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cause a portion of the wing lift to be lost. One way to compensate for this 
effect is to analyze the configuration at 1/2 degree greater incidence than that 
of the experiment. This was done for the 0141 and 05 comparisons shown 
in the results section. 


For some specialty applications or for cases where computing costs are not 
a concern, excellent results have been obtained by using 150 crude iterations 
followed by 150 cycles in the crude/fine system. Other investigators have 
improved flow simulations by increasing the relaxation factor from w = 1.5 to 
u> = 1.75. The new user should begin by using standard cycle/factor values 
and experiment upward if convergence or comparisons are not satisfactory. 


5 . WING TIP DI VERGENC E 

Analyses performed on highly loaded wings with severe taper and high 
sweep indicate that occasionally a diverging numerical condition can develop at 
the wing tip. It occurs only when very strong shock waves form at the wing 
tip and experience indicates that occurrences are close to conditions for which 
the flow separates. The numerical condition is caused by improper zones of 
dependence and/or unfavorable wing tip mesh cell aspect ratios. A rotated 
difference scheme is now under development. The new scheme eliminates this 
problem but requires additional refinements and an extensive check-out effort. 
A simpler "fix" has been achieved by automatically reducing wing tip mesh 
resolution. This resolution reduction scheme has been automated by computing 
a "probability of divergence" factor which is based on Mach number, inci- 
dence, wing thickness, and planform taper. Occasionally, an application might 
surface for which the automatic system fails. For these cases, the user can 
specify the number of fine mesh points at the wing tip (the input variable 
AQXFG) . The highest resolution possible without divergence indicators 
(pressure "wiggles", increasing should be sought. 


6. COMPARISONS IN HIGH SPANWISE GRADIENT REGIONS 

At transonic speeds, flow gradients might be high in the wing spanwise 
direction. This is particularly noticeable near wing-fuselage junctures, 
wing-winglet junctures, or mid- wing regions where multiple shock wave sys- 
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terns coalesce. The computational method predicts chordwise pressure dis- 
tributions at a number of spanwise locations. These "output" stations may not 
be aligned with experimental data stations or locations where output is 
required. Since small span position discrepancies might result in significant 
pressure differences, some applications may require that output data be inter- 
polated or extrapolated. 

On a related matter, it should be noted that the fuselage side is modeled 
by mesh points positioned along a wing span cut. Depending on the wing/ 
fuselage spacial relationship, the "physical" body side might be as much as 1/2 
mesh cell away from the "computational" body side. This presents no problem 
for most applications, but for some cases where exact juncture pressures are 
required, and the 1/2 mesh cell shift occurs, the flow simulations may be 
compromised. 

7. GRID GENERATION 

The computational method has been developed with the philosophy that 
engineering users should not have to modify grid systems for most applica- 
tions. Instead, the systems should be constructed in the same consistent 
manner for each case. It can be seen by studying the input data format that 
only the configuration geometry and flow conditions are input. Grid gen- 
eration has been automated. A knowledgeable user might modify the coordinate 
systems for special applications, but this requires sufficient familiarity with 
the coding to make changes to the FORTRAN. 

This automated grid philosophy does not guarantee that an "optimum" 
system will always be generated, particularly for complex configurations with 
many components. For example, a transport might feature two nacelles posi- 
tioned below each wing at different heights. The code in its present form is 
capable of representing only a single height or position for both nacelles. 
Thus it may be necessary to perform two different analyses to obtain results 
for each component at its proper position. 

An attempt has been made to insure that a nacelle or pod is modeled at 
its proper distance from the wing by adjusting the grid so that grid points 
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representing the nacelle centerline are positioned at the requested nacelle 
location. This does not insure, however, that the computational nacelle wing 
gap will equal the physical gap. This is an important aspect of many flow 
simulations. The user should observe what the computational gap is between 
the points which represent the side of the nacelle facing the wing and the 
points which represent the wing plane. If this gap is not equal to the phys- 
ical space between the components, the nacelle height variable should be 
adjusted until a match is obtained. 

8. SUPERCR ITICAL AIRFOIL/ WING COVE SEP A RAT ION 

The boundary layer method incorporated is of the finite difference type 
which terminates once separation is encountered. The viscous analysis is 
initiated after an initial inviscid pressure field is generated. For super 
critical- type sections the initial severe unfavorable pressure gradient found in 
the cove or lower surface trailing edge region tends to induce a premature 
separation. A temporary displacement thickness ( 6 *) must be "fit" for the 
cove region to allow the solution process to continue. This fitting process is 
quite crude and often a true, computed, attached flow boundary layer is never 
achieved despite the fact that experimental data indicates the flow is attached. 

It has also been noted that the code tends to predict supercritical section 
cove pressures that are higher than experimental values. Predicted wing 
section lift levels might be several percent too high with an attendant nose 
down pitching moment increment. 

9. SURFACE SHAPE IRREGULARITIES 

The present method does not employ conformal mappings for wing sections 
and no smoothing operation is performed on wing section ordinates which are 
input. No attempt should be made to "smooth" surface breaks caused by 
control surface deflections or irregularities resulting from poor manufacturing 
processes. Experience indicates that the code will simulate the resulting 
pressure disturbances quite well. For some applications, these predictions will 
be useful. 
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10. WINGS WITH STRONG VISCOUS INTERACTIONS 

Some applications might feature wing contours for which viscous effects 
are quite severe. Good agreement can be obtained for these cases, but the 
total number of iterations may have to be increased by a factor of 3 or 4. 


143 


144 


I 



REFERENCES 


1. Murman, E. M., and Cole, J. D. : "Calculation of Plane Steady Transonic 

Flow," AIM J. , vol. 9, Jan. 1971, pp. 114-121. 

2. Bailey, F. R. ; and Steger, J. L. : "Relaxation Techniques for Three- 

Dimensional Flow about Wings," AIM J. , vol. 11, Mar. 1973, pp. 318-325. 

3. Keller, J. D. ; and South, J. C. , Jr.: "RAXBOD: A Fortran Program for 

Inviscid Transonic Flow over Axisynmetric Bodies," NASA TM X-72831, 
Feb. 1976. 

4. Jameson, Antony; and Caughey, D. A.: "A Finite Volume Method for 

Transonic Potential Flew Calculations." A Collection of Technical 
Papers - AIM 3rd Computational Fluid Dynamics Conference, June 1977, 
pp. 35-54. (Available as AIM Paper 77-635.) 

5. Boppe, C. W. ; and Stern, M. A.: "Simulated Transonic Flows for Aircraft 

with Nacelles, Pylons, and Winglets," AIM Paper 80-0130, Jan. 1980. 

6. Boppe, C. W. : "Transonic Flow Field Analysis for Wing-Fuselage 

Configurations," NASA CR-3243 , May 1980. 

7. Hesse, W. J. ; and Mumford, N. V. S. : Jet Propulsion for Aerospace 

Applications , Pitman Corp. , NY, 1964, p. 128. 

8. Liepmann, H. W. ; and Roshko, A.: Elements of Gas Dynamics , John Wiley and 

Sons, Inc., NY, 1967, p. 53. 

9. Freuler, R. J. ; and Gregorek, G. M. : "An Evaluation of Four Single 

Element Airfoil Analytic Methods," Advanced Technology Airfoil Research 
Conference, NASA CP-2045, vol. I, part 1, Mar. 1978, pp. 133-162. 

10. deKuyper, R. E. : "Wind Tunnel Tests of a 1/8.8 Scale Reflection Plane 

Model of the Grumman Gulf stream III Airplane," CA1SPAN Report No. M-4001- 
W-26, Nov. 1977. 


PAGE PLANK NOT FILMS!) 


145 



11 . 


Flechner, S. G. ; and Patterson, J. C. , Jr.: "Tabulated Pressure Measure- 
ments cm a Large Subsonic Transport Model Airplane with High-Bypass-Ratio, 
Powered, Fan-Jet Engines," NASA Till X-2530, May 1972. 


12. Montoya, L. C. ; Flechner, S. G. ; and Jacobs, P. F. : "Effect of 

on a First-Generation Jet Transport Wing, II - Pressure and 
Load Distributions for a Semispan Model at High Subsonic Speeds 
TN D-8474, July 1977. 


Winglets 
Spanwise 
." NASA 




1. Repoft No. 2. Government Accession No. 

NASA CR-4066 


4. Title and Subtitle 

Aerodynamic Analysis for Aircraft With 

Nacelles, Pylons, and Winglets at 

Transonic Speeds 


7. Authors} 

Charles W. Boppe 


9. Performing Organization Name and Address 

Grumman Aerospace Corporation 
Bethpage, NY 11714 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


3. Recipient's Catalog No. 


5. Report Date 
April 1987 


6. Performing Organization Code 


8. Performing Organization Report No. 

10. Work Unit No. 

11. Contract or Grant No. 

NASl-14732 

13. Type of Report and Period Covered 

Contractor Report 

14. Sponsoring Agency Code 

534-02-13-13 


15. Supplementary Notes 

Langley Technical Monitor: Perry A. Newman 

Final Report , Part I 

16. Abstract 

A computational method has been developed to provide an analysis 
for complex realistic aircraft configurations at transonic speeds. 
Wing-fuselage configurations with various combinations of pods, 
pylons, nacelles, and winglets can be analyzed along with simpler 
shapes such as airfoils, isolated wings, and isolated bodies. The 
flexibility required for the treatment of such diverse geometries 
is obtained by using a multiple nested grid approach in the finite- 
difference relaxation scheme. Aircraft components (and their 
grid systems) can be added or removed as required. As a result, 
the computational method can be used in the same manner as a wind 
tunnel to study high-speed aerodynamic interference effects. 

The multiple grid approach also provides high boundary point I 

density /cost ratio. High resolution pressure distributions can J 

be obtained. Computed results are correlated with wind tunnel I 

and flight data using four different transport configurations. 
Experimental/computational component interference effects are I 

included for cases where data are available. The computer code • 

used for these comparisons is described in the appendices. 


17. Key Words {Suggested by Author(s)) 18. Distribution Statement 

Transonic Flow 

Computational Aerodynamics Unclassified — Unlimited 

Flow Simulations 


Finite Difference Relaxation 
Interference Effects 

19. Security Classif. (of this report) 20. Security Classif. {of this page) 

Unclassified Unclassified 


Subject Category 

21 No. of Pages 22. Price 

154 A08 


02 


it" 


I 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA^Langley, 198/ 













